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ABSTRACT 

The star formation rates (SFR) of low-metallicity galaxies depend sensitively on the gas metallicity, 
because metals are crucial to mediating the transition from intermediate-temperature atomic gas to 
cold molecular gas, a necessary precursor to star formation. We study the impact of this effect on the 
star formation history of galaxies. We incorporate metallicity-dependent star formation and metal 
enrichment in a simple model that follows the evolution of a halo main progenitor. Our model shows 
that including the effect of metallicity leads to suppression of star formation at redshift z > 2 in dark 
halos with masses < 10 11 M , with the suppression becoming near total for halos below ~ iq 9 - 5-10 
Mq. We find that at high redshift, till z ~ 2, the SFR cannot catch up with the gas inflow rate 
(IR), because the SFR is limited by the free-fall time, and because it is suppressed further by a lack 
of metals in small halos. As a result, in each galaxy the SFR is growing in time faster than the IR, 
and the integrated cosmic SFR density is rising with time. The suppressed in situ SFR at high z 
makes the growth of stellar mass dominated by ex situ SFR, meaning stars formed in lower mass 
progenitor galaxies and then accreted, which implies that the specific SFR (sSFR) remains constant 
with time. The intensely accreted gas at high z is accumulating as an atomic gas reservoir. This 
provides additional fuel for star formation in 10 10 — 10 12 M & halos at z ~ 1 — 3, which allows the SFR 
to exceed the instantaneous IR, and may enable an even higher outflow rate. At z < 1, following the 
natural decline in IR with time due to the universal expansion, the SFR and sSFR are expected to 
drop. We specify the expected dependence of sSFR and metallicity on stellar mass and redshift. At a 
given z, and below a critical mass, these relations are predicted to be flat and rising respectively. Our 
model predictions qualitatively match some of the puzzling features in the observed star formation 
history. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: high-redshift — galaxies: 
ISM — ISM: molecules — stars: formation 



1. INTRODUCTION 

This paper attempts a first-principles theoretical anal- 
ysis of metallicity-dependent star formation as a ma- 
jor player in determining the cosmological star-formation 
history. To help motivate this study, it is useful to start 
with a review of the observational context and the ten- 
sion between these observations and current models. 



Th e ACDM cosmological model (Blumenthal et al. 
19841 has been very successful in explaining the large- 
scale structure of the universe, and the general prop- 
erties of dark-matter halos in which galaxies form, in- 
cluding their mass distribution as a function of time and 
environment. While we can crudely connect the halos 
to observed ob j ects via statistical inferen ce (e.g. Con- 
roy et al.l|2006| |Conroy fe Wechsler|[2009| , we stiirTack" 
a first-principles model for determining the observable 
stellar content of a given dark-matter halo. Numerical 
simulations and semi-analytic model that attempt to do 
so generally encounter a number of problems. 

First, the cosmic star formation rate (SFR) density is 
observed to be rising in time from z > 8 till z ~ 3, and 
it remains high till z ~ 1 (e.g. Hopkins fc Beacom||2006 
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Kistl er et al.|2009| |Bouwens et al.|2010[ |Horiuchi fe Bea- 



com|2010p , whileACDM predicts a strong decline of cos 
mological in flow rate (IR) into ha los of a given mass at all 
epochs (e.g. |Neistein et aT||2006[ ). Since the star forma- 
tion timescale is short compared to the Hubble time for 
most of cosmic history, the SFR and the IR are expected 
to be tightly locked - galaxies convert gas into stars at a 
rate that is ultimately determined by the gas supply rate. 
As a result, in most current models, the decline in IR at 
high z drives a decline in SFR density that starts much 
earlier than z ~ 1 — 2 (e.g.|Murali et al.||2002| |Hernquist 
I fc Springel|2003l [Springel fc Hernquist|2003||Keres et al 
2005 1 |Baugh et al.|2 005; Bo uche et al.|2010|ISchaye et al 
|2010[ ). There is a related observational indication for a 
rapid growth in SFR within a given population of halos as 
they grow in time, fast er than the predicted mass growth 
by accretion till z ~ 2 (Papovich et al. 2010 ). The tension 
between observation and theory is emphasized by the ob- 
servational indications for a "sSFR plateau", where the 
specific SFR (i.e. the SFR per unit stellar mass) is rather 
constant between z ~ 8 and z ~ 1 — 2, and is significantly 
lower than predicted at z > 6 (Weinmann et al. 2011 



and references therein). The above may also be relatec 
to the "missing dwarf problem" , where the low observed 
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abundance of small galaxies compared to the model pre- 
dictions indicates that the models tend to form stars too 



rapidly at high rcdshift ( 


Fontanot et al. 2009 


Marchesini 


et al. 2009 Cirasuolo et aT2010 


) . All the above indicate 



SFR at low masses and at high redshift. 

At intermediate redshifts, z ~ 1 — 2, the models seem 
to face an opposite problem. The observed SFR and 
sSFR are high, and perhaps eve n exceed the predicted 



insta ntaneous IR and specifi c IR ( Daddi et al. 2007 Dave 
2011 Weinmann et al.|2011). The problem is made even 



worse by observational indic ations for massive outflows 
from galaxi es at z ~ 2 — 3 (Steidel et al. 2010 Genzel 



et al.||2011[ ), which in some cases may exceed both the 
SFR and the predicted instantaneous IR, thus indicating 
that the freshly accreting gas is not the only source of 
gas supply. 

The need to suppress SFR in small halos and at high 
redshift has been addressed in models by the introduction 
of various "feedback" mechanisms, which are assumed 
to inhibit star formation and/or remove gas from galax- 
ies. Proposed mechanisms include p hotoionization by 
the UV background (see the review by Loeb fc B arkanaj 
2001[ ), and stellar feedback by supernovae ( i Dekel &; STTkl 
1986[) or radiation pressure (Thompson et al |20 05 



Mur- 



ray et al.||2010| |Hopkins et al.| |20lff The attempts 
made so tar to implement stellar feed back in simulations 
and s emi-analytic models (SAM) are not really sat i sfac- 
tory ([Springe! h Hernquist||2003[ |Kobayashi||2004[ |Cen | 
&; Ustriker||200t>| |Uppenheimer fe L>ave||200tj| |(Jovernato| 



et al.||2007p , because they operate on scales too small 
to be resolved, and have to be incorporated via simplis- 
tic recipes that sometimes fail to capture the complex 
sub-grid physics involved. As a result, the modeling in 
most cases is limited to an attempt to tune the feedback 
prescriptions in order to obtain a better fit to the ob- 
servation. Since different sub-grid prescriptions produce 
different results, they have limited predictive power. 

The feedback mechanisms mentioned above are as- 
sumed to be effective in removing the gas from halos 
of virial velocity V v < 50 kins -1 . This is consistent 
with the indications for lack of gas in such systems at 
z r~j 2 — 3, based on observations of damped Lyma n a ab- 



sorption systems o n the lines of sight to quasars ( Barnes 



& Haehnelt 2010, and references therein). However, m 
order to match the obs erved SFR history a nd other ob 



servational constraints, Bouche et al. ( 2010[ ) have argued 
that the SFR has to be suppressed also m more massive 
halos, in the mass range M v ~ 10 10 — 10 n M@, corre- 
sponding at z ~ 2.5 to V v = 50 — 100 kms -1 . Halos in 
this mass range do seem to contain g as at the level of half 



the c osmological baryonic fraction (Barnes & Haehnelt 



2010 1, which calls for a non-ejective feedback mechanism 
The observational indications for SFR and outflow rate 
that exceed the IR at z ~ 1 — 3 also call for a feedback 
mechanism that causes the accumulation of gas reservoirs 
in smaller galaxies at higher redshifts. 

In this paper we consider from first principles a phys- 
ical process that may offer an entirely different solution 
to the problems of the star formation history. In the past 
decade, observations of star formation in nearby galaxies 
have unambiguously revealed that star formation is as- 
sociated with the molecular phase of a galaxy's interstel- 
lar medium (ISM) ( |Wong & Blitz|2002| [Kennicutt et ah 



2007| |Leroy et al.|2008| |Bigiel et al.|2Q08 ) . In gas where, 
averaged over kpc scales, the dominant chemical state is 
H I rather than H2, the star formation rate is more than 
an order of magnitude lo wer than in molec ular gas of 
comparable surface density (jBigiel et al. 2010 1. It is likely 



that even in such H i-dominated regions star formation is 
only occurring in a transiently-formed molecular phase. 
Beyond the local universe, observations at high redshift 
indicate a similar phenomenon: H i-dominated damped 
Lyman-a systems and outer parts of Lyman break galax- 
ies both show star formation rates more than an order of 
magnitude lower than would be ex pected for molecular 



gas of comparable surface densities ([W olf e fc Chen|2006 
Wild et al |[2007} | fiafelski et al.||2010) . 



Krumholz et al.| ( |2011 ) show that this observation can 
be understood in term s of the therm odynamics of in- 
terstellar gas (also see Schaye 2004 1. The H 1 to H 2 



transition occurs in gas at high volume and column den- 
sity; a high volume density produces a large H2 formation 
rate, while a high column density produces greater shield- 
ing against dissociating interstellar UV photons. The 
gas temperature depends on volume and column density 
in nearly the same way, because high volume densities 
raise the rate of cooling by collisionally-excited metal 
lines, while high column densities increase the amount of 
shielding against grain photoelectric heating. As a result, 
over many orders of magnitude in gas volume density, 
column density, and metallicity, the H 1 to H2 transition 
is an excellent proxy for the transition from gas that is 
too warm to form stars (T > 100 K) and gas that is cold 
enough (T ~ 10 K) to undergo runaway collapse to stel- 
lar densities^] The su bsequent numerical simulations of 
Glover & Clark (2011) confirm this picture. 

The distinction between the warmer atomic and colder 
molecular phases of the ISM is significant even in the 
present-day universe, but it is crucial in the early uni- 
verse. The transition between the phases is highly sensi- 
tive to the metallicity of the gas, occurring at a very dif- 
ferent column density in the Milky Way, at Solar metal- 
licity, than in t he Small Magellanic Cloud, at ~ 20% So- 
lar metallicity (|Tumlinson et al.|2002||Bolatto et al.|2007l 
Leroy et al.||2007[ |Krumholz et al.||2009b[ |Bolatto et al j 
2011J) . In low metallicity blue compact dwarf and dwarf 



irregular galaxies, which are thought to be local analogs 
of common high-redshift star-forming galaxies, H 1 sur- 
face densities reach ~ 100 Mr pc~ 2 (Fumagalli et al. 



2010 1, whereas in Solar metallicity systems a transition 
to H2 generally prevents H 1 surface dens ities from ex- 



ceeding 10-20 M Q pc~ 2 (Bigiel et al. 2008 ). Since metal 



licities are lower in the early universe, the metallicity- 
dependence of the warm H 1 to cold H 2 transition, and 
thus of the star formation rate, has a potentially dra- 
matic effect on cosmic star formation history. 

Despite these observations, only a few simulations and 
models of star formation over cosmological time have in- 

1 A note on terminology: in many cosmological applications it 
is common to refer to gas at temperatures ~ 10 4 K as cold to dis- 
tinguish it from hotter ~ 10 6 K halo or IGM gas. Our terminology 
here will be closer to that used in the ISM literature, where cold 
refers to gas at typical molecular cloud temperatures of ~ 10 K, 
and warm refers to typical H I temperatures of ~ 10 2 — 10 3 K. 
The distinction is important because it is not possible to produce 
Jeans masses as low as ~ Mq unless the gas can cool to ~ 10 K 
temperatures. 
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corpo rated the effects of metal licity-dependence. Analyt- 
ically, Krumholz et al. ( 2009a ) show that observations of 
the column density and metaliicity distribution of DLAs, 
and their lack of star formation, can be explained by 
modeling the w arm H I to cool H 2 transition . On the 



simulation s id e, [Ro bertson fc Kravtsov (2008), Gnedin 
et al.| (I2009[ ), |Gnedin & Kravtsov| ( |2010j ), and |Kuhlen 
et al. ( 201 ip present simulations showing that low metal- 
iicity systems, either dwarfs in the local universe or star- 
forming galaxies at high redshift, form stars at a much 
lower rate than found in earlier simulations omitting 
metaliicity effects. However, thus far these simulations 
have only been able to sample relatively small cosmolog- 
ical volumes, and have not been able to advance past the 
peak of the star formation history at z ~ 2. Thus they 
convey limited information about the large-scale effects 
of metalli city-dependent s tar f ormation. Semi-an alytic 
models by |Fu et al!] ( |2010[ ) and |Lagos et al. ( 2010[ ) have 
begun to explore how "metallicity-aware" star formation 
laws affect galaxy populations^ However, these models 
have numerous moving parts, and it is unclear exactly 
how the results depend on the numerous other parame- 
ters already built into the SAM. 

In this work we present a simple toy model to under- 
stand how the need to undergo a warm H I to cold H2 
phase transition as a prelude to star formation might 
affect the cosmic star formation history. Due to the nu- 
merous approximations we make, we do not expect this 
result to produce quantitatively exact results. Our goal 
is simply to understand the basic physics of the effect 
and its qualitative implications. The model can be sum- 
marized in three basic evolution equations for a single 
halo, which we present in Section [2l We then describe 
the results of this model in Section |3J and we compare 
these results to observations in Section 01 We summarize 
our conclusions and discuss implications in Section [5] 

2. HALO EVOLUTION MODEL 

In this section we present a simple model for the evo- 
lution of the main progenitors of dark matter halos and 
their baryonic content. We characterize a halo by its to- 
tal mass Mh and by the gas, stellar, and metal masses 
M g , M*, and Mz in the disk galaxy that is assumed to 
form at the halo center. The mass of metals here includes 
only metals in the gas, not metals locked in stars. 



2.1. Go 



id Stellar Mass Evolution 



Foll owing the "bath-tub" model of Bouche et al 
(2010), we assume that the gas and stellar mass in a 
halo are determined by three processes: inflow, star for- 
mation, and ejection of material by star formation feed- 
back. The continuity equations that govern this system 
are 



M„ = M n 



(1 - R)M t 



form 



M* = M» >ln + (l-#)M*, form . 



(1) 
(2) 



Here M gjin and M, 



arc the rates of cosmological inflow 

of gas and stars (S ection 2.1.1 1, Af*,f rm is the star forma- 
tion rate (Section |2.1.2 ), H is the stellar return fraction, 



2 |Obreschkow et al.| ||2009| and |Obreschkow fc Rawlings| |2009| ) 
also studied the hi I to M2 transition in their semi-analytic models, 
but they included it only in post-processing, so the star formation 
rate was not affected by the molecular fraction. 



and e ou t is the amount of mass ejected into the IGM per 
unit mass of stars formed (Appendix [A]) . We summarize 
the parameters that appear in these equations, as well 
as others that appear below, in Table [TJ Below we dis- 
cuss our estimates for the quantities appearing in these 
equations. 

2.1.1. Accretion rates 

The average total (baryonic plus dark matter) inflow 
rate into a halo of virial mass Mh at redshift z due to 
cosmological instrcaming from the cosmic web is approx- 
imated by 

M Ma = -aM^dj, (3) 

where M M2 = M h /10 12 M , a = 0.628, f3 = 0.14. Here 
uj = 1.68/ D(t) is the self-similar time variable of the EPS 
formalism, in which D{t) is the linear fluctuation growing 
mode, and it is approximated by 



-0.0476[1 



0.093(1 



-1.2212.5 



Gyr- 1 . (4) 



This a ppr oximation has been deri ved by Neistei n et al.| 
( |2006[ ) and |Neistein fc DekelM2008| ) based on EPS theory 
and was fine-tuned using the Millennium cosmological N- 
body simulation. We adopt the cosmological parameters 
from WMAP: fl m = 0.27, Q A = 0.73, h = 0.7, and 
o~ g = 0.8l(^] This expression is accurate to better than 
5% for z = 0.2 - 5, and to - 10% for z = - 10. The 
mass dependence is accurate to 5% for halos from mass 
Mh t i2 = 0.1 — 10 2 . The distribution of accretion rates 
about this mean can be fit by a lognormal with a stan- 
dard deviation of ~ 0.2, plus a tail representing m ergers 



that extend s to a factor of 10 above the average (Dekel 
et al.|[2009 Goerdt et al., 2011, in preparation), though 
we do not implement this scatter in our model. Although 
we use equation ^ for all our numerical computations, 
a simpler expression, which we will use to make scaling 
arguments later, is 



M h 



34M yr- 1 M, 1 ; 1 1 4(i + z )2 



(5) 



This matches equation Q to ~ 10%. We note that if 
the power 1.14 is approximated by unity, and if the 2.4 is 
replaced by 5/2 (which is the exact predicted value in the 
Einstein-deSitter cosmology phase that is approximately 
valid at z > 1), this accretion rate corresponds to a mass 
growth rate Mh cx exp (—73), with 7 ~ 0.9. 

The corresponding gas and stellar inflow rates into a 
halo are 



M„ 



- 0.17ei n /j(, fi.17 fg.ixiMh 



M* )i n = 0.17/ M .i 7 (l-/ 9> i n )M/ 1 , 



(6) 
(7) 



where fh is the cosmic baryon fraction, fb,a.n — 
fa, in/0. 17, / g ,in is the fraction of the inflow that reaches 



J The values of a and /3, and the numerical coefficients in equa- 
tion Bl given here have been updated to th is cosmology, and are 
therefore slightly different than those given in |Neistein etal] |2006 ) 
or |Neistein fc Dekel| ( |2008^ . Also note that formally equation H in- 
cludes all accretion, both smooth and in the form of major mergers. 
One might therefore worry that the median accretion rate might 
be less than the mean given by the equation. However, since major 
mergers are strongly subdominant for the ranges of halo mass and 
redshift we be will considering, this is a relatively small effect. A 
more quantitative discuss ion of the expected leve l of variation in 
accretion rates is given in |Neistein fc Dekel ( 2008]) . 
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the halo as gas rather than stars, and e;„ is the frac- 
tion of accreted gas that reaches the galactic disk rather 
than being shock-heated and going into the galactic halo. 
Note that we have not attempted to include the effects 
of a reduction in accretion rates onto small halos with 
virial velocities due to suppression of cooling by the UV 
background after the universe is reionized. This is ex- 
pected to be a large effect for halos with virial velocities 
below 30 km s" 1 , and a factor of a few effect for halos 
with virial velocities of 30 — 50 km s _1 (Thoul & Wein- 
berg 1996). Our omission of the effect will make some 
difference for the very smallest halos that we model, but 
the majority of the models we present below are in the 
regime where suppression of accretion is unimportant or 
is only a factor of a few effect. 

In our fiducial model we compute / Sl i n using a self- 
consistent approximation described in Appendix [51 and 
we explore how this affects our results by also considering 
a variant model in which we take / Sl m = 1 for all halos 
at all redshifts. 

We have estimated the fraction of gas inflow from out- 
side the virial radius that actually reaches the galaxy at 
the halo center, ej n , using analytic arguments and hydro 
cosmological simu l ations (|Birnboim fc Dekel|2003| |Keres| 
et al.| 2005i|2009b||Dekel k ijirnboim|^006| |Ocyirk et al.| 
Deke l et al.||2009p . Well below a, critical halo mass 
1 there is no virial shock, and the gas flows in 
so e; n ~ 1. At z > 2, in halos above the critical 
mass, cold streams bring in most of the gas along the 
dark-matter filaments of the cosmic web, and they pene- 
trate effici ently through the hot medium, yielding e m ~ 1 
also there ( Dekel et al.|20 09). A statistical analysis of 400 
AMR-simulated halos at z ~ 2.5 indeed reveals e- m ~ 0.9, 
with the exact value depending on how the averaging 
is performed (in pre paration) . Another estimate based 



2008 
cold, 



on SPH simulations ( Faucher-Giguere et al. 2011) gives 
similar values at high z but somewhat smaller values, 
e; n ~ 0.5, at z ~ 2.5. At z < 2, in halos more massive 
than the shock-heating scale, the cold flows are broader 
and their penetration power is weaker, so the accretion 
rate into the disk is suppressed ( [Dekel fe Birnboim|2006 



Ocvirk et al. 2008). This turns out to cxplain the evolu- 



tion of the red sequence o f galaxies ( Dekel fc Birnboim 
2006 Cattaneo et al.|2006 ). We crudely model this using 



a simple analytic form fol 



owing 



Dekel et al. (2009): 



Mi2,max = max 2, 



1 



£in = 



3M* )PS (z)y 

1, M M2 < M 12 .max 
0, M M2 > Mia .max- 



(8) 
(9) 



Here M* ) ps(z) is the Press-Schechter mass at a given 
redshi ft, which we compute from the |Sheth fc Tormen] 
(1999) ellip soidal collapse mode l following the procedure 
outlined in Mo & White (2002). We have also run mod- 
els in whi ch we adopt "the analytic fitti ng formulae pro- 
posed by Faucher-Giguere et al. (2011), and found that 
the results are qualitatively unchanged. Of course the 
use of a sharp mass cutoff for accretion is an oversim- 
plification of reality. Clearly some galaxies with masses 
below Mi2 jmax are early types without accretion or star 
formation, while some with larger masses are late types 
that have ongoing star formation. Indeed, based on cos- 



mological simulations, the transition from a population 
of galaxies dominated by full cosmological accretion rate 
to a population dominated by shutdown is stretche d over 
more than a decade in halo mass (e . g. Figure 1 of [ Birn- 
boim et al.|2007| |Ocvirk et aL|2008[ |Keres et al.l|2009a[ ). 
However, given the simplicity of the remainder ot our 
model, there would be little point in attempting to in- 
clude this scatter. Moreover, given the difficulty in car- 
rying high resolution simulations of cosmological inflow 
to redshifts z%l, theoretical estimates of the accretion 
rates there must be considered highly uncertain in any 
event. In this paper we will be more concerned with 

redshfits 2^2, where equation |9| may be regarded as 
reasonably accurate. 

2.1.2. Metalhcity-Dependent Star Formation 
The mo del presented in the p revious section is similar 



to that of |Bouche et al. (2010). We now diverge from 
that model by taking into account the atomic-molecular 
phase transition in computing the rate at which gas turns 
into stars. On scales ranging from individual clouds to 
entire galaxies, the star formation rate is well-described 
by ~ ef[MH 2 /iff, where Mjj 2 is the total molecular 
gas mass, tg is the free-fall time of the gas, and eg ~ 0.0 1 
dKrumholz fc Tan||2007[ |Krumholz fc Thompson] ! 2007 ) . 
.'his result is an integrated version of the classical Ken- 



nicutt 



( 1998 ) relation based on observations 
However galaxies at low surface densities and metallic 
ities are observed to form stars at a rate considera bly be 
low that pred i cted by the Kennicutt relation (e.g. 
fc Che n||2006l IWild et al.|[2007| |Leroy et al.j 2008 
308 TW" ... 



Wolfe 
Bigic 



et al.||20"08[ |Wyder et al.||2009p , a result that has~TTeen 
successfully explained by noting that stars form only in 
molecular gas, and that systems with low surface density 
and metall icity tend to have little of their gas in molec 
ular form ( Robertson fc Kravtsov [20081 iKrumholz et al. 
2009a|. ' 



Gncdin et al. 2009; Gncdin & Kravtsov 2010 



We therefore adopt a star formation rate m a given halo 



M, 



* .form 



= 2tt 



(10) 



where S g is the gas surface density, and /h 2 is the frac- 
tion of that gas in molecular form. All quantities inside 
the integral are functions of the distance from the disk 
center. 

For the surface density, we assume that the gas in disks 
follows an exponential profile S = 'E c e~ r / Rd . The scale 
length Rd is assumed to be proportional to the virial 
radius of the halo in which it resides, 



R d = 0.05A .ii?„ 



(11) 



where A = (l/v / 2)( J/M h )/(VR) is the spin parameter 
for a halo of angular momentum J, mass Mh, virial ra - 
dius R v , and circular velocity V ( |Bullock et al. 2001), 
and Ao.i = A/0.1. Note that we choose to make Rd half 
of XR V because for an exponential disk the half-mass ra- 
dius is 1.7 scale radii, and given the uncertainties our 
ansatz places the half-mass radius at roughly XR V . Tidal 
torque theory and N-body simulations give on av erage 
A ~ 0.04 but the observed radii of z ~ 2 galaxies (Gen 
zel et al. [20061 12008[ ) suggest that the g. 
somewhat higher angular momentum. I 



ms in them 



has a 



Following Dutton 
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et al. (2010), we adopt A ~ 0.07 as our typical value. 
1'he virial radius is related to the virial mass and the 
expan sion factor a = 1/(1 
2006 1 

Rv.ioo 



z) by (Dekel & Birnboim 



1.03M^ 3 2 A 1/3 , 



(12) 



where i?„,i o = Rv/100 kpc, A 1/3 = .4/(1/3), 

— 1 /3 

A = (A 20 0^m,0.3^0.7) °> ( 13 ) 

and fi m .o.3 = ^m/0.3. Here 

A(a) ~ (18tt 2 - 82Q A {a) - 39fi A (a) 2 )/Q m (a) (14) 

is the overdensity within the virial radius at a given 
epoch, A 20 o = A/200, fi TO (a) = Q m a- 3 /(Q A + tt m a~~ 3 ), 
and Q\(a) — l~ft m (a). At redshifts > 1, in the Einstein- 
de Sitter regime, A200 ~ 1 and so A ~ a for our standard 
cosmology. In this cosmology A 2 oo grows to 1.7 at red- 
shift 0. Given this surface density profile, the central 
surface density is related to the disk gas mass by 



£ C ,< 



2nRj 



0A25M gtll \^M-% 3 A-? 3 , 



(15) 



where £ c ,o = £ c /(l g cm 2 ). Note that, if M g oc Mh, 
the surface density scales with mass and redshift like oc 

M 1 /3( 1 + z )2_ 

Also note that our treatment of disk sizes is extremely 
simple. We have ignored effects like the variation of halo 
concentration with mass and redshift, and the possible 
evolution of galactic angular momenta with reds hift or 
other galaxy properties (e.g. Burkert et al. 2010). The 
disk size affects our models mainly by changing the sur- 
face density and thus the radius at which gas transitions 
from atomic to molecular, as discussed in the next sec- 
tion. This transition is also affected by metallicity, and 
the uncertainties associated with metallicity evolution 
are probably at the order of magnitude level. Given the 
size of these uncertainties, a more detailed treatment of 
disk radii seems unnecessary. 

The local free-fall time d epends on the gas surface den- 
sity £ of the galactic disk. Krumholz et al. ( 2009c I have 
analyzed the structure of molecular clouds and developed 
a theoretical model for star formation that gives 



f « (2.6 Gyr)" 1 
iff 



(Eo/0.18)" - 33 , 
(£ /0.18) - 33 , 



£ < 0.18 
£ > 0.18 ' 

(16) 

where £0 = gem 2 ). The low £ regime cor- 

responds to galaxies like the present-day Milky Way 
where molecular clouds are confined by self-gravity, while 
the high £ regime describes galaxies like low-redshift 
ULIRGs where they are confined by external pressure. 
This formula agrees well with the star formation rate 
observed in nearby g alaxies (c.f. Figures 1 and 2 of 
Krumholz et al.|[2009c ), and we adopt it here. 

Most importantly, the fraction of the ISM in the cold, 
molecular phase, /h 2 , is determined by the balance be- 
tween grain photoelectric heating and UV photodissoci- 
ation on one hand, and collisionally-excited metal line 
cooling and H2 formation on dust grains on the other 
hand. (Near zero metallicity other processes become im- 
portant, but, as we discuss below, we will not consider 
this regime.) |Krumholz et al.| ( |2008"} |2009b| ) and |McKee" 



& Krumholz (20101 show that these processes produce a 
molecular traction that depends primarily on the surface 
density and metallicity of the galactic disk, and relatively 
little on any other parameters. A crude approximation 
to their result is that the molecular fraction is 



/h 2 ~ £/(£ + IQZ^Mq p<f 



(17) 



where £ is the total gas column density and Z$ = 
(Mz/M g )/Z@ is the metallicity normalized to the So- 
lar neighborhood value Z Q = 0.02. Thus regions with 
£ <C 10 Zq 1 Mq pc -2 are primarily atomic, and those 
with £ S> WZq 1 Mq pc~ 2 are primarily molecular. A 
more accurate expression, which we adopt here, is 



/h 2 



1 



\l+0. 25s J ' 



S < 2 

v 0, s>2 
ln(l + 0.6x + 0.01x 2 ) 



X = 3.1- 



0.6r c 

1 + z - 365 



4.1 

t c = 320cZ £ , 



(18) 

(19) 

(20) 
(21) 

where c is a clumping factor that accounts for smoothing 
of the surface density on scales larger than that of a single 
atomic-molecular complex. If £ is measured on 100 pc 
scales then there is no averaging and c » 1 , while on ~ 1 
kpc scales it is ~ 5 ( Krumholz et al. 2009c I , and we adopt 
this as a fiducial parameter in our models. This analytic 
model agrees very well with numerical simulations that 
follow the full ti me-dependent chemistry of H 2 formation 



and dissociation (Krumholz & Gnedin 



2011). 



It is also possible to form stars in galaxies that are es- 
sentially free of molecules and dust, a topic that has been 
studied exte nsive ly since the pionee ring work of |Bromm| 
et al.| ( |2002[ ) and [Abel et al.| ( |2002[ ). In these cases the 
cooling required tor star formation is driven by the tiny 
fraction of H 2 that is able to form by gas phase processes, 
rather than via grain catalysis, the dominant process over 
most of cosmic time. However, this population III pro- 
cess is extremely slow and inefficient compared to the 
normal mode of star formation. While it is important for 
providing the seed metals that enable the normal mode 
to begin, we will neglect the contribution of population 
III stars to the total star formation rate of the Universe. 

2.2. Metallicity Evolution 

The final ingredient to our model is the metallicity of 
the gas, which affects its ability to form a cold molecular 
phase. The first seed metals that begin the process of H 2 
formation and the transition to normal star formation 
are produced by these population III stars. A single pair- 
instability supernova from one of these stars pollut es its 
host halo u p to a metallicity of ~ 10~ 3 of solar (Wise 
We therefore adopt Zjqm — 2 x 10 _o as 
"intergalactic medium" metallicity at which 
all halos start, and which characterizes newly accreted 
gas. We test our sensitivity to this choice below by also 
considering a model with Zigm ten times higher. 

Once the molecule fraction becomes non-zero and nor- 
mal star formation begins, we set the metallicity based 
on the expected yield of the resulting stars. Following 
the standard practice in chemical evolution models (e.g. 



et al. 


2010) 
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Rcdshift evolution of the total halo mass (solid blue lines), total baryonic mass M g + M* (dashed red lines), gas mass 



Fig. 1. 
M g 

Mh 12 = lj an d M-hVX = 10, as indicated at the top of each panel. All masses are plotted in units of Mi2 
we also show the halo mass multiplied by the universal baryon fraction /j, (dashed blue lines). 



0.008, M h ,%2 = 0.016, 
M/10 12 M Q . For comparison 



Maeder 1992), we use the instantaneous recycling ap- 



proximation and write the metal production rate in terms 
of the yield y of the stellar population. The net metal 
production rate corresponding to a star formation rate 
M*,form is 2/(1 — R)M* t { otm ; we adopt y = 0.069, as dis- 
cussed in Appendix [A] However, only a fraction of those 
metals will be retained in the galaxy rather than lost in 
supernova explosions. This fraction is a very sharply in- 
creasing function of halo mass, and based on a rough fit 
to the simulations of Mac Low & Ferrara ( 1999 ) we take 
the fraction of metals ejected by supernovae to be 



C = CloC 



2/M re 



(22) 



where M rot is the halo mass (in units of 10 12 M Q ) at 
which supernovae become unable to eject most of their 
metals, and (j is the fraction of metals retained in halos 
much smaller than M rpt . As a fiducial value, based on 



the simulations of Mac Low & Ferrara (19991), we adopt 
Mret — 0.3. Choosing a fiducial value of Q is more diffi- 
cult, because even if almost all metals are init ially ejected 
from the galactic disk in low mass halos (as Mac Low & 
Ferrara find) , many of these will not escape the halo en- 
tirely, and will later be re-accreted - as must be the case, 
since halos below Mh,\2 = 0.3 are not completely devoid 
of metals. In the absence of firm theoretical guidance we 
adopt as a fiducial value (j = 0.9, i.e. we assume that 
90% of metals will be ejected completely from small ha- 
los, while 10% will be retained or re-accreted. Below, we 
test the sensitivity of our results to both M ret and (i by 
varying these fiducial values. 

The two other factors capable of changing the metal 
content are accretion of pristine material, which adds 
metals at a rate Z\Qy[M g - ln , and expulsion of existing 
interstellar medium by feedback, which removes metals 
at a rate (Mz / Mg)e ou tM*f orm . This expression assumes 
that the expelled gas has a metallicity equal to the mean 
metallicity of the galaxy. Combining metal production, 
accretion of pristine material, and ejection of existing 
interstellar material, the evolution equation for the mass 
of metals is 



• = tf(l-fl)(l-C)-M,,fenn 



M z ■ 

+ ZiGMMg i n — e ou t—— M* form- (23) 
Ma 



Note that we distinguish between the ejection of hot met- 
als that are already part of a 10,000 km s _1 SN shock and 
the driving of outflow in the cold dense ISM, the former 
being easier to eject and more strongly dependent on the 



halo mass via £. We have verified that our results are 
rather insensitive to the value of e ou t and to a possible 
halo-mass dependence in it; for more discussion of this 
issue, see Appendix |A"| 

Finally, we caution that our treatment of metals ne- 
glects the existence of metallicity gradients within galac- 
tic disks. In the local universe, these are observed to be 
relatively modest. For spiral galaxies the center-to-edge 
metallicity difference is typically at most half a dex (e.g. 
Oey & Kennicutt 1993 ), with the metal licity gradient dis - 
appearing completely outside R25 (e.g. Werk et al.|2011 ). 
Dwarf g alaxies also have ess entially no metallicity gradl- 
ent (e.g. Croxall et al.| [2009). Thus metallicity gradients 
seem unlikely to very significantly alter our results. Any 
potential effects are certainly likely to be smaller than 
those induced by changing ^igm by an order of magni- 
tude, as we do below. 

3. MODEL RESULTS 

With the model ingredients described in Section [2] in 
place, we can now calculate the evolution of star forma- 
tion in growing main progenitor halos. The total mass, 
gas mass, and stellar mass then evolve according to equa- 
tions ([3]) , , and ( 2]) , while the metal mass evolves fol- 
lowing eq uati on ( [23 ) . The star formation rate is given by 
equation (|10|, and the gase ous infall fraction / Sl i n is cal- 



culated from equation (B4). In addition to the fiducial 



model, for comparison we repeat also compute several 
variants of the model. All model parameters are summa- 
rized in Table [TJ 

In the first variant, we set /h 2 = 1 everywhere regard- 
less of column density or metallicity. This allows us to 
isolate the effects of the requirement that H2 form before 
stars do on the star formation history of the universe. In 
the second, we set f g ,i n = 1, so that the inflow contains 
no stars. This allows us to understand the importance of 
these processes. In the third variant, we use M ret = 0.03 
instead of 0.3 as the characteristic mass at which halos 
start retaining most of their metals. In the fourth variant 
we set Cio = 0.8, thereby allowing halo that are below the 
threshold to retain 20% of their metals rather than 10% 
as in the fiducial model. Variants three and four allow us 
to test how changing our recipe for the ability of halos to 

retain metals, eq. ( 22 1, alters our results. In the fifth vari- 

-4 



ant we change the IGM metallicity to Zigm = 2 x 10" 
which is 10 times the fiducial value, and corresponds to 
1% of Solar. The true metal content of the IGM is likely 
between these ex tremes and probably evolves with red- 
shift as well (e.g. Schaye et al. 2003 Simcoe|[2011| , but 
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Fig. 2. — Redshift evolution of the star formation rate M t f orm (solid black lines), gas inflow rate M g \ n (dashed red lines), and stellar 
inflow rate M*,i a (dashed black lines) for the same halos as in Figure^ 
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Fig. 3. — Redshift evolution of the specific star formation rate sSFR = M, f orm /M* (solid black lines), specific gas inflow rate sIR = 

M t ,i n /M t (dashed black lines) for the same halos as in Figure [l] 

that stars formed in mergers account for ~ 10% of the 



^s,in/^» (dashed red lines), and specific stellar inflow rate sIR* 



for simplicity we consider only these two models, which 
should allow us to bracket reality. Finally, in the sixth 
variant, which we refer to as the Egp = 10 M pc~ 2 
model, we adopt 



/h 2 



E < 10M©pc- 2 
E > lOAf pc~ 2 



(24) 



i.e. where we allow star formation to occur only where 
the surface density exceeds 10 M© pc -2 independent of 
metallicity. This is very similar to the treatment of star 
formation adopted in many numerical simulations and 
semi-analytic models that do not include the physics of 
the atomic-molecular phase transition. 

For both our fiducial model and its variants, we 
compute a set of 400 model halos starting at z = 
30 and ending at z = 0. The halos have initial 
masses uniformly spaced in \ogMh,i2 from logM/^12 = 
—5.54 to log Mh, 12 = —4.06, corresponding to present- 
day masses \ogMh,\2 = —3 to 3. A halo of ini- 
tial mass M^ init begins its evolution with gas mass 
M ff)in it = Q.17/6,o.i7/o,init- M /»,init, stellar mass Af*,i„ it = 
0.17/(,, .17(1 - / ff ,init)Mh,mit, and metal mass M z = 
ZiGuMg^Tut, where we set / 9i i n it = 1 for the variant 



model with f g 



1, and / fl ,init = 0.9 for all other mod- 



els. The choice / g ,i n it deserves some comment. It is 
clearly not realistic to assume that halos begin entirely 
devoid of stars (/ g ,init = 1), an d that they form them 
only in the quiescent mode that our calculation includes. 
In mergers, or in the accretion shocks that appear when 
halos first begin to accrete baryons, gas surface densities 
can reach values much higher than our exponential disk 
model would predict, and these high surface densities 
will allow some star formation even in halos with very 
low metallicities. However, it is not entirely clear what 
value of / g ,init we should choose to represent this effect. 
Our choice of 0.9 is motivated by observations indicating 



cosmic star formation rate budget at high z ( Rodighiero 
et al.|20lT I . While this contribution is obviously not the 
same for all halo masses, as a crude estimate we sim- 
ply assume that all halos convert at least 10% of their 
baryons to stars early in their lives. 

In Section 3.1 we present results for a f ew se lected halos 
using our fiducial model, and in Section [3. 2| we compare 
these to our variant models. We present the population 
statistics of our halos in Section 1331 



3.1. Fiducial Model 

We consider four example halos, with present-day 
masses of M M2 = 0.008, 0.016, 1.0, and 10.0; these 
masses are chosen to illustrate several interesting behav- 
iors in the model. We note that each of these halos has 
a virial velocity above 25 km s _1 at all redshifts below 
z <~ 10, and above 30 km s _1 from z = 2 — 3, when the 
bulk of the gas is accreted and most star formation oc- 
curs. Thus these halos will not be significantly affected 
by by ionization-driven winds ( Shaviv fc Dekel||2003[ ), 
though for the smallest ones the accretion rates may be 
reduced by a factors of a few due to reduction in cooling 
by the UV background (|Efstathiou|1992| |Thoul fc Wein- 
berg|1996[[Hoeft et al.|2006l ). Mote, however, that below 
we will plot results for smaller halos for which photoion- 
ization effects are likely be important. As noted above, 
we have not attempted to include these simply to avoid 
further complicating the model.) Figures [l] - [4] summa- 
rize the evolutionary history of our four example halos 
using the fiducial model. Figure [l] shows the mass of each 
component, Figures [2] and 3] show the total and specific 
rates of star formation (including only quiescent star for- 
mation, not stars formed in accretion-driven bursts), as 
well as gas and stellar inflow, and Figure [4] shows the 
evolution of the gas-phase metallicity and column den- 
sity. 
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Fig. 4. — Redshift evolution of the metallicity normalized to Solar ZJZq = (Mz /M g ) /Zq (solid black lines) and central column density 
S c = M g / (2-kR"^) (dashed black lines) for the same halos as in Figurelll For comparison, we also show the column density S mo i at which 
a gas with the given metallicity would become 50% molecular (dasheared lines), as computed from equation (TTSp. 



3.1.1. Star formation cannot catch up at high redshift 

We see that, at very high redshift (z ^> 5) all halos are 
dominated by gas rather than stellar mass, and they are 
forming stars more slowly than they are accreting both 
gas and stars from the intergalactic medium. However, 
there is star formation during this phase, which is critical 
for building up metals in the ISM. Although the halos' 
starting mctallicitics are quite low, this is compensated 
for by their high column densities, and as a result much 
of their mass is able to form molecules and participate in 
star formation. They remain gas-dominated simply be- 
cause their star formation rate is unable to keep up with 
the inflow rate. To see this, note that the characteristic 
timescale for a halo to double its mass via inflow is (using 
the approximate equation [5] for the mass inflow rate) 



= 2.1 Gyr Af££(l + z) 3 



-2.4 



(25) 



where (l+z) 3 = (l+z)/3. The timescale for gas inflow to 
change the gas mass is similar during the era when star 
formation has not yet significantly reduced the gas mass. 
In comparison, for a halo with M g = 0.17 'fb,o.n fgMh, 
where f g is the fraction of the baryonic mass in gas, we 
can estimate the time required to turn order unity of 
the molecular gas in a galaxy into star s by plugging M g 
into equation (1 1 5 h and then equation (16). This gives a 



timescale for star formation at the disk center 
eg 



tsF 



= 2.5 Gyr A°;f fr»;™f s 



0.33 t-0. 33j^-0.lri , „\-0.67 



g ""h,12 (1t*73 ' 

(26) 

where we have taken A fa a and we have assumed that 
we are in the S Ci o > 0.18 regime; both of these approxi- 
mations hold well at high z. Thus we see that, for z 2, 
iacc <C tsF, so we expect that the star formation rate will 
not be able to keep up with the inflow rate, more so at 
higher redhifts. This effect is even further enhanced by 
the fact that only a portion of the gas is cold and molec- 
ular. We note that even if igp scaled with the disk and 
halo crossing times, which scale with the Hubble time 
and thus oc (1 + z) -15 , the SFR was still unable to catch 
up with the inflow rate at high enough redshift. We 
should emphasize that the inability of the SFR to catch 
up with the I R is a result of having the disk radius given 
by equation (11). During a major merger the gas can 



be driven to a much smaller radius, leading to a shorter 
free-fall time and more rapid star formation. However, 
there are not enough major mergers to make a substan- 
tial difference in gas consumption on cosm ological scales 
( |Neistein fc Dekel||2008| |Dekel et al.||2009[ ). 



3.1.2. Metal buildup and the peak in SFR 

Despite the inability of the SFR to catch up with the 
IR, however, star formation does take place and build 
up metals. More massive halos have larger surface densi- 
ties, and as a result their ratio of star formation rate and 
metal production rate to new gas inflow rate is higher. 
Furthermore, the y retain a larger fraction of their met- 
als (Equation [22]) . Both effects lead more massive halos 
build up metals more quickly than less massive ones dur- 
ing this phase, as seen in Fig. [4] 

As time passes halo surface densities drop due to the 
rise in disk scale lengths that accompanies expansion of 
the universe, and this has divergent effects on the differ- 
ent halos. The lowest mass halo, with a present-day mass 
= 0.008, only builds up metals very slowly, and so 
the surface density at which it undergoes the atomic to 
molecular transition is relatively high. At z ~ 8, its cen- 
tral surface density is high enough to produce a ~ 50% 
molecule fraction, but by z ~ 6 it falls too low to produce 
many molecules, and the star formation rate suddenly 
drops. This becomes self-reinforcing, since new inflow 
continues to bring in pristine gas, and without a con- 
tinuous source of metal production within the halo the 
gas-phase metallicity drops back toward the IGM value. 
The halo effectively shuts down. 

In the halo with a present-day mass Mh.12 — 0.016, 
on the other hand, metal production is more rapid, and 
the galaxy remains able to form molecules and stars in 
its center for its entire evolutionary history. Moreover, at 
z ~ 1 — 3, the gap between the central surface density and 
the surface density required to have a 50% molecular frac- 
tion increases (see Figure [4]). As a result, an increasing 
fraction of the disk becomes star- forming r] which allows 
the star formation rate to rise significantly faster than 
the inflow rate, until it slightly exceeds the inflow rate 
near z ~ 2. This burst of star formation at a rate above 
the gas inflow rate contributes to the peak in the star 
formation history of the universe, as we will see below. 

In contrast, the two highest mass halos have higher 
metallicities and column densities. As they produce more 
and more metals, the ratio of their central surface densi- 
ties to the atomic-molecular transition density rises. This 
allows molecules to form and star formation to take place 
over more and more of their area as we approach z ~ 2. 



4 Qualitatively, such inside-out growth is a robust prediction 
of our model that is consistent with the star formation histories 
inferred for nearby g alaxies using resolved stellar populations (e.g. 
|WilIiams et al.|2009| , but a quantitative comparison will require a 
more detailed model than the simple one we have computed. 
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While their star formation rates are still nearly an order 
of magnitude below their accretion rates at z ~ 5, by 
z ~ 2 they have produced enough metals for their star 
formation rate to match their accretion rate. This effect 
helps suppress star formation at high z and enhance the 
peak at z ~ 2. 

Perhaps the most striking conclusion so far, as can be 
read from Figure[2j is that our model predicts for the red- 
shift range z = 8— 2 a steep increase with time of the SFR 
within each growing massive galaxy. It could be approx- 
imated by M*,form oc exp(— 0.65z). This is very different 
from the assumption sometimes adopted by modelers of a 
SFR that is decaying exponentially with time (the "tau" 
model). In this redshift range, the accretion rate itself 
is growing slowly with time, reflecting the exponential 
growth of halo mass in equation [5] The SFR is growing 
much more steeply, partly because iaccAsF < 1 at high 
z and is increasing with time, and partly because the 
metallicity-dependent quenching is strong at high red- 
shift and is weakening with time (see below). By z ~ 2, 
the SFR catches up with the accretion rate as the latter 
has slowed down and the metallicity has become suffi- 
ciently high. As a result of the quenching of SFR at 
high redshift, as seen in Figure |T gas keeps accumulat- 
ing at z > 4, and it provides additional supply for star 
formation at 2^4—2. 

3.1.3. Drop beyond the peak 

Past z ~ 2, the accretion and star formation rates in 
both the Mh,i2 — 1 and Mh,\2 = 10 halos drop. This 
drop is a result of the competition between growth in 
the accretion rate as halo mass increases at fixed redshift 
(Mh oc M^ ) an d decline in accretion rate with redshift 
at fixed halo mass (Mh oc (1 + z) 2A ) as a result of the 
expansion of the universe. At high z halos grow fast 
enough for the former effect to dominate, so accretion 
rates rises as one approaches the present. Past z ~ 2, 
however, the decline in accretion rate associated with 
expansion of the universe begins to dominate, and the 
accretion rate falls. For the Mh,i2 = 10 halo this decline 
is further accelerated when the halo reaches a mass of 
~ 2 and gas accretion shuts off. By z = the star 
formation rate in this halo has dropped off from its peak 
value to nearly zero. In contrast, the Mh,i2 = 1 halo has 
only declined in star formation rate by a factor of ~ 10. 

3.2. Comparison to Variant Models 
3.2.1. The effect of metallicity on SFR history 

In Figure [5] we compare the SFR in our fiducial model 
to the SFRs in our variant models. We see that the vari- 
ants and the fiducial model are all very similar for the two 
higher mass halos. These are able to form molecules with 
relatively little difficulty due to their consistently high 
column densities, so their star formation histories are not 
significantly affected by metallicity-dependent star for- 
mation. The story is different for the lower mass halos. 
For them making the accretion purely gaseous (/ 9 ,i n = 1) 
or reducing the metal retention mass (M ret = 0.03) still 
have relatively modest effects - the star formation lasts 
slightly longer in the Mh.12 = 0.008 halo, and peaks 
slightly sooner in the Mh.12 = 0.016 one - but quali- 
tatively there is little difference. In contrast, assuming 
that gas can form stars regardless of whether it is in a 



cold, molecular phase (/h 2 = 1), that this phase tran- 
sition occurs at a metallicity-independent threshold of 
Esf = 10 Mq pc~ 2 , decreasing the maximum fraction of 
metals from supernovae that are ejected (£io = 0.8), or 
increasing the IGM metallicity (Zigm = 2 x 10~ 4 ) have 
dramatic effects in these small halos. All these changes 
allow much more rapid star formation at high 2, reduce 
the star formation peak behavior at z ~ 1 — 2 in the 
Mh,i2 = 0.016 halo, and allow star formation to continue 
up to the present day in the M/, 12 = 0.008 halo. Thus 
we see that including a realistic treatment of ISM chem- 
istry and thermodynamics in models is likely to have a 
dramatic effect on the star formation history of the uni- 
verse, and to make star formation sensitive to exactly 
how metals are distributed in halos, while changing the 
stellar inflow fraction or the rate of star formation in 
accretion shocks has relatively little effect. 

The insensitivity of our results to the choice of M let 
might at first seem somewhat surprising, given how dra- 
matically metallicity-dependent star formation changes 
the halos' SF history. However, note that the halos where 
metallicity makes a difference have masses much smaller 
than M r e t even if we reduce M ret by a factor of 10 com- 
pared to Mac Low fc Ferrara s value. This suggests that 
what matters in determining the mass below which SF 
will be suppressed is not the precise mass at which ha- 
los start retaining most of their metals, but instead the 
ability of halos below this retention mass to build up a 
non-negligible level of metallicity, which in turn depends 
on the ratio of metal production rate to accretion rate 
in these halos. This is a function of the star formation 
rate, the yield, the metallicity of the accreting gas, and 
the fraction of metals that are retained in small halos0 

3.2.2. The origin of a sSFR plateau 

In Figure[6j we show the specific star formation rate for 
our fiducial model and for all the variants. Here we see 
that setting f g ^ n = 1 has a major effect on the results. It 
is easiest to understand this by focusing on the two higher 
mass halos, since for these both the fiducial model and 
the variants have about the same total star formation 
rate. We find that the fiducial model and all the variant 
models except / g ,i n = 1 have specific star formation rates 
in the two higher mass halos that are relatively flat at 
high z. This is fairly easy to understand. At redshifts 
above z ~ 4—5, in all of these models the stellar accretion 
rate exceeds the in situ star formation rate (c.f. Figures 
[2] and [3). By examining where in Figure u\ the stellar 
mass curve begins to deviate from simply tracking the 
total baryonic mass, one can see that stars formed in 
situ do not dominate the total stellar mass until z ~ 3. 
Thus, at 2 > 3, the stellar mass in a halo of mass M^ is 
approximately M* = 0.17/(,,o.i7(l — fg,in)Mh- Assuming 
the gas mass is M g — Q.YJ fb,o.nfg,\uMh (i.e. that gas has 
not been significantly depleted by star formation yet) , we 
therefore have a specific star formation rate 



sSFR 



AL 



forn 



/H: 



(27) 



M* 1 — fg t i n V"" tff ( 

5 Although we have not directly investigated how varying the 
yield might change our results, note that the yield y enters the 
evolution equations only through the combination y(l — £). Thus 
changing the yield should be nearly equivalent to changing (j , and 
the yield is uncertain only at the tens of percent level, as compared 
to factors of order unity for the factor 1 — f . 
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Fig. 5. — Redshift evolution of the star formation rate M. M f orm in the same halos as in Figure 1 for our fiducial model {solid black lines), 
for our variant models: fn 2 = 1 (dashed black lines), f g ,i n = 1 (dashed red lines), M re t = 0.03 (dot-dashed green lines) (j = 0.8 (dashed 
blue lines), Ziqm = 2 x 10 — 4 (dot-dashed purple lines), and Egp = 10 Mq pc~ 2 (dashed aqua lines) . 




Fig. 6. — Redshift evolution of the specific star formation rate 
the same models as in Figure [5] 

where the angle brackets indicate a surface density- 
weighted average over the disk. Since the latter is gener- 
ally ~ (2 — 3 Gyr) -1 , we obtain a nearly constant specific 
star formation rate of 1 — 3 Gyr -1 at high redshift. The 
emergence of a sSFR plateau at high redshifts is insen- 
sitive to the exact value of f g ,im &s long as the fraction 
of ex situ stars is non- negligible, but it affects the am- 
plitude of the sSFR plateau, and its extent toward lower 
redshifts. Metallicity-dependent star formation flattens 
the sSFR as a function of z even further, because it makes 
high z halos even more dominated by accreted / accre- 
tion shock- formed stars than in the model with /h 2 = 1. 
In contrast, in the models with f g ; m = 1 halos are dom- 
inated by in situ stars at all epochs and masses. As a 
result, is very small at high z, giving rise to a larger 
specific star formation rate. Since tsF gets close to t acc 
as we approach the present day, M* rises faster than Mh , 
and the specific star formation rate declines toward the 
present, rather than remaining flat. 

At intermediate redshifts, z ~ 2 — 4, where the in-situ 
SFR is already a major contributor to the stellar mass, 
the fact that the SFR is growing roughly exponentially 
with time leads to a similar exponential growth of its 
integral, the stellar mass, and thus helps extending the 
plateau toward z ~ 2. 



sSFR = M,. form/-^* i n the same halos as in Figure^ Lines represent 

M* fiducial 



3.3. Population Statistics 

We now turn our attention to the ensemble properties 
of our entire model grid. Figure [7] shows the star forma- 
tion rates and gas accretion rates in the halos as they 
evolve in mass with redshift. To help interpret this plot, 
in Figure [8] we show two ratios: 



fiducial 



(28) 



sir; 



M, 



(29) 



g, in, fiducial 



The first of these describes the factor by which the SFR is 
reduced by our inclusion of metallicity quenching relative 
to a model that does not include. The second describes 
the factor by which the SFR is suppressed relative to the 
rate of gas inflow in our fiducial model. 

These plots show that a star formation law in which 
stars form only in cold molecular gas dramatically sup- 
presses star formation in low mass halos. Notice that 
contours of constant SFR and large sh 2 are close to hor- 
izontal at z > 2, indicating that metallicity quenching 
acts much like a halo mass threshold for star formation 
at z > 2. The value of this threshold mass is sensitive 
to how well small halos are able to retain their metals 
and to the metallicity of accreting IGM gas, changing 
from ~ 10 10 M & in our fiducial model to ~ 10 9 Mq in 
our variant where small halos retain twice as much of 
their metal production, and to slightly smaller masses in 
the model where the IGM metallicity is raised to 1% of 
Solar. As we will see below, this induces a ~ 0.5 dex 
shift in the total star formation rate budget of the uni- 
verse at high z. In contrast, both the gas inflow rate and 
the star formation in the model where stars form equally 
well in warm atomic or cold molecular gas (/h 2 = 1) 
or where the atomic-molecular transition is taken to be 
metallicity- independent (Ssf = 10 Mq pc~ 2 ) continue 
all the way down to the smallest halos in our model set. 
For halos above the threshold, the story is quite different. 
At z > 4, metallicity quenching reduces the SFR in these 
halos too, but only by factors of sh 2 ~ 2. In contrast, the 
SFR lags the gas inflow rate by factors of sir ^3 — 10, 
simply because star formation cannot keep up with the 
inflow rate. Conversely, at z < 4 we find that sh 2 < 1 
in massive halos, indicating that star formation rates are 



11 













— -0.5 -| 







SFR, f„ 



SFR. M„ 



0.03 



1.0 



SFR, 



c.s 



-;■ 



O 



SFR, E, 



10 



-2.0 



-3.0 



2 4 6 8 10 

z 



Fig. 7. — Star formation rate (SFR) and gas inflow rate (IR) as a function of halo mass and redshift. The rates are indicated by color, 
and are limited to our model grid. Heavy white or black lines are contours of constant SFR or IR, running from log(M / Mq yr~ x ) = —1.5 
to 2.5 in steps of 1 dex, as indicated. Light gray lines indicate the mass versus redshift for 20 sample halos in our grid of 400. The panels 
show the SFR in our fiducial model and in our variant models and the IR in our fiducial model, as indicated. 



actually higher in the fiducial model than in the /h 2 = 1 
or Esf = 10 Mq pc~ 2 models. The reason for this be- 
havior is that metallicity quenching reduces the SFR at 
high z when the metallicity is low, but it does not remove 
the gas. At lower z when the metallicity rises, this gas is 
able to form stars. Thus metallicity quenching in large 
halos simply delays star formation, helping to produce 
the peak of cosmic SFR at z ~ 2. In contrast, the choice 



.in 

or M r , 



0.03 clearly makes little difference. We 



therefore conclude that the star formation history of the 
universe is sensitive to the need to form a cold phase, 
and that the ability of different halos to do so depends 
on their ability to retain metals. The SF history of the 
universe does not depend strongly on stars being formed 
externally or in accretion shocks or on the exact mass at 
which halos become efficient at retaining metals. 

Figure [9] shows the corresponding specific star forma- 
tion and specific inflow rates. We see that the fiducial 
model exhibits a very wide range of halos masses and 
redshifts over which the specific star formation rate is 



~ 2 Gyr -1 . At lower z the specific star formation rate 
starts to fall, as star formation both increases the stellar 
mass and decreases the gas mass, but this process takes 
some time to complete, and in most halos the specific 
star formation rate does not drop below 1 Gyr -1 un- 
til after z = 2, and does not fall below 0.5 Gyr -1 until 
z *~ 1. The model with /h 2 = 1 exhibits a slightly larger 
range of specific star formation rates, but except for very 
low mass halos, it is qualitatively similar to the fiducial 
model. Thus we conclude that molecules do not dramat- 
ically change the specific star formation rate in massive 
halos. 

In contrast, the model with f g ^ n — 1 reaches sSFRs of 
many tens per Gyr, and the specific inflow rate reaches 
similarly large values. This is the same effect we iden- 
tified in the previous section in studying the histories of 
individual halos. When f g> \ n =/= 1, halos at very high z 
tend to be dominated by stars that are directly accreted 
or formed in accretion shocks, so their stellar mass sim- 
ply scales with the halo mass. Since the star formation 
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Fig. 8. — SFR suppression factor in our fiducial model relative to the model without metallicity quenching (left) an d re lative to the 
inflow rate (right). The format is similar to Figure[7] but here colors indicate the factor sh 2 (equation |28[ l or sir (equation [29jl. The former 
is the factor by which metallicity quenching suppresses star formation, and the latter is the ratio by which the SFR is smaller than the 
gas inflow rate in our fiducial model. Contours indicate values of logs = —1.5, — 1, —0.5, 0, and 0.5 as indicated, i.e. suppression factors 
of 0.032, 0.1, 0.32, 1, and 3.2. Note that s > 1 corresponds to an enhancement of the star formation rate rather than a suppression. Also 
note that in the upper left of the plot (high M^, low z) the SFR and IR in all the models become very small, and as the result the value 
of s is set by numerical noise in the integrator. To suppress this effect in this plot we set s = wherever the fiducial model SFR is below 
10~ 3 Ms 



rate scales with the gas mass (and thus also with the halo 
mass), the specific star formation rate is constant. For 
fg.in = 1, on the other hand the stellar mass is no longer 
proportional to the halo mass at high z, and the sSFR 
becomes large. 

It is interesting to notice the contrasting roles played 
by metallicity quenching (i.e. the fact that /h 2 7^ 1) and 
stellar inflow (i.e. the fact that f g ^ n 7^ 1). In order to 
have a SFR that increases with time but a sSFR that 
is flat, as the observations discussed in Section [I] appear 
to demand, one requires that the SFR deviate from the 
inflow rate, but that the total stellar mass does follow 
the time-integrated inflow rate. These seemingly contra- 
dictory requirements are met by a combination of sup- 
pression of star formation by metallicity quenching and 
increase in stellar mass due to stellar inflow. 

4. COMPARISON TO OBSERVATIONS 

In this section we use our model grid of 400 halos to 
generate predictions for a number of observed correla- 
tions between various properties of galaxies that have 
been reported in the literature, or that will become ob- 
servable in the next few years as new facilities come on- 
line. Unless stated otherwise, all the comparisons in this 
section use our fiducial model. 

4.1. SFR in a growing galaxy 

We can also use our models to predict the evolution of 
SFR or other galaxy properties with observational sam- 
ples selected in a variety of ways. One observational ap- 
proach is to select galaxies based on their comoving num- 
ber density. In this strategy, at each redshift z one selects 
galaxies with luminosities L near a cutoff luminosity L(z) 
chosen based on the condition that the comoving number 
density of galaxies with luminosities L > L{z) satisfy the 
condition that n{> L{z)) = no, where uq is some fixed 
number density. This approach is useful because the co- 
moving number density of main halo progenitors should, 
to good approximation, remain constant with redshift. 
Thus this technique should produce a selection that is 
close to following one of our theoretical main progeni- 
tor tracks in Figure [7j It is therefore a straightforward 



prediction of our model. 

To compare our models to samples of this sort, wc 
first compute the halo abundance n(M^,z), the num- 
ber density of halos at redshift z w ith masses in the 
range M h to M h + dM h , using the |Sheth fc Tormen] 
(1999) approxima tion, following the co mputation proce- 
dure outlined by |Mo fc White ( |2002| ). At each z, we 
then numerically solve tor the halo mass Mh{z) for which 
/mh(z) n (Mh, z) dMh = no, where no is the target num- 
ber density. We then interpolate between the two model 
halos in our grid that most closely bracket Mft,(z)r] 
In Figures 10 and 11 we compare the SFR ana stel- 



lar mass as a function of lookback time for our model 
halos to three number densit y- selected samples , taken 



from Marchesini et al. 



Papovich et al. (2011) 



([2009 > 

and Lundgren et al. (2012, in preparation). The observa- 
tional samples use a threshold density no.obs = 2 x 10~ 4 
Mpc -3 for Papovich et al. and Marchesini et al., and 
1.8 x 10 -4 Mpc -3 for Lundgren et al. For the latter 
sample, only stellar masses are available, not star forma- 
tion rates. The choice of no for the theoretical compari- 
son requires some care, because depending on the selec- 
tion method the sample is likely to be at least somewhat 
incomplete. Halos at a given mass will show some scat- 
ter in their luminosity, and downward scatter (e.g. a low 
UV luminosity associated with a temporarily low SFR) 
will remove galaxies from the sample. Upward scatter 
can also add lower mass galaxies to the sample, but the 
effects are not symmetric because, if star formation his- 
tories are bursty, then galaxies will tend to spend more 
time with SFRs below their l ong-term average than with 



SFRs above their average, van Dokkum et al. (20081 



estimate that ~ 50% of galaxies at z ~ 2 will have sup 
pressed SFRs, so that the appropri ate value of np f or a 
UV-selected sample such as that of |Papovich et al7| is a 

6 In deriving our comparison from the observational grid in 
this manner, we implicitly assume that luminosity or stellar mass 
(which is used to derive the number density in the observed sam- 
ple) is a monotonically increasing function of halo mass. This as- 
sumption is satisfied by our model grid, at least over the redshift 
ranges where luminosity or stellar mass selection are used for the 
observations. 
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Fig. 9. — Specific star formation rate and gas inflow rate as a function of halo mass and redshift. This is the same as Figure 
here colors indicate specific rather than total star formation rate or inflow rates. The heavy white and black contours indicate speci 
formation or gas inflow rates of sSFR or sIR = 0.05, 0.5, 1, 2, and 4 Gyr , as indicated. 

4.2. Evolution of sSFR at Fixed Stellar Mass 



but 

star 



factor of ~ 2 above the nominal value. The samples of 
Marchesini et al. are stellar mass selected and those of 
Lundgren et al. are H-band selected, and should there- 
fore be substantially more complete; for these samples 
the appropriate value of no to choice is probably close to 
the nominal one for the survey. Thus we plot curves for 
our model galaxies selected using no, model = 2 x 10~ 4 , 
4 x 10~ 4 , and 8 x 10~ 4 Mpc~ 3 , corresponding to per- 
fect completeness, factor of 2 incompleteness, and factor 
of 4 incompleteness in the observations. As the figures 
show, for reasonable estimates of observational complete- 
ness, the models do a good job reproducing the observed 
trends in SFR and stellar mass. In particular, we see 
that in both the models and the data, SFRs and stel- 
lar masses rise roughly exponentially with —z at high z 
before leveling off near z ss 2 — 3. 

In fact, the evolution of SFR and stellar mass in the 
model halo with Mi 2,0 = 10 shown in Figures [JJ and [2j 
which has a halo mass of 10 11 Mq at z ~ 6, is very similar 
to the observed evolution seen in Figures [10] and [TTJ 



Another common observational method is to examine 
the properties of gala xies at fixed stellar mass across a 
range in redshift (e.g. i Papov ich et al.||2006| |Stark et al 
20101. To 



2009] [Gonzalezet 
tions of this sort 
10 9 



al.||2010) . ' To compare to observa- 
we pick a target stellar mass in the 
10 11 M m for which high-redshift obser- 



range ~ iu — ±v ivi & 
vations are available, and at each z in our model grid 
we select the halo with stellar mass closest to the target 
value. We then plot the corresponding sSFR in Figure 
"12] For comparison, we also plot a region illustrating the 
range of observationally-de termined sSFR values b ased 
on a data compilation from Weinmann et _al. (2011 ). We 
see the same general result as in Figure [9j in our fidu- 
cial model the sSFR has a nearly constant value of ~ 2 
Gyr -1 from z ~ 8 to z ~ 2. This is in good agreement 
with the observations. In contrast, the specific gas in- 



fall rate sIR = M gM /M„ 
star formation rate if / g ,i n 
magnitude higher at z ~ 



in these halos, and the specific 
1, is more than an order of 
than at z ~ 2, in disagree- 
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Fig. 10. — Star formation rate in density-selected galaxies as 
a function of lookback time, normalized to the lookback time at 
redshift z = 11. This approximates the evolution of SFR in a 
grow ing main-progenitor g alaxy. Squares are observations taken 
from |Papovich et al.||2011[ ) selected to have observed number den- 
sity no = 2 X 10 — * Mpc -3 at all redshifts. Error bars are their 
suggested values. Lines are values taken from our fiducial model, 
selected to have number densities no = 2 X 10~ 4 , 4 X 10~ 4 , and 
8 X 10 -4 Mpc -3 , from top to bottom. The smallest value would 
be appropriate if all the central ga laxies were forming stars at a 
high rate and the Papovich et al.| sample is perfectly complete, 
while the other two values assume that the sampled star-forming 
galaxies are 50% or 25% of the central galaxies. For details on how 
we derive the model values, see the main text. 
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Fig. 12. — Specific star formation rate sSFR = M*,form/-M* in 
our fiducial model (solid lines) and in the f g - m = 1 model (dotted 
lines), and specific gas accretion rate sIR = M g ; n /Af* (dashed 
lines) at fixed stellar masses of M* g = 1, 5, and 20 (black, red, 
and blue lines, respectively). The M, 9 = 20 line does not extend 
past 2 ~ 6 because even the highest mass of our grid of model halos 
has not reached this stellar mass at higher redshifts. We do not 
show the /jj 2 = 1 or Af re t = 0.03 models because they does not 
differ significantly from the fiducial case in this halo mass range. 
For comparison, we also show the observed sSFR (shaded region). 
The r egion shown is a fit to the data compiled by (Weinmann ct al. 
|2011| see their Figure 1). 
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Fig. 11. — Same as Figure |10| but we now plot the stellar 
mass rather than the star formation rate. Squares again represent 
UV-selected galaxies. C ircles are stellar mass-selected data from 
|Marchesini et ah] | |2009[ |, and triangles are H-band selected data 
(which should be similar to stellar mass-selection) from Lundgren 
et al. (2012, in preparation). The error bars on the UV and stellar- 
mass-s elected data sets are the values suggested by Papovich ct al. 
j2011| . Data points at very similar lookback times in fact repre- 
sent samples at the same lookback time, but have been offset from 
one another for clarity. 

ment with the observations. The underlying reason for 
this trend is the same as for the individual galaxies: at 
high z, stellar mass is roughly proportional to halo mass 
because most stars are accreted, and star formation rate 
is also roughly proportional to halo mass because isF has 
only a weak redshift dependence. This provides an expla- 
nation for the otherwise puzzling observed sSFR plateau, 
as discussed in ^T]and |5] 

4.3. The Star- Formation Sequence 

Yet another important observable is the correlation be- 
tween stellar mass and star formation rate for star form- 




FlG. 13. — The star formation sequence, stellar mass versus 
specific star formation rate, at 2 = (green), 2 = 1 (lavender), 
2 = 2 (red), 2 = 3.7 (blue), 2 = 6 (indigo), and 2 = 8 (orange), 
from bottom to top, in our model grid (solid lines). Lines end at the 
maximum stellar mass sampled by our model grid at that redshift. 
For co mparison, hatched regions show t he observed correla tions at 
2 = princhmann et al.||2004[ l, 2 = 1 plbaz et al.|2007) , 2 = 2 
| |Daddi er"aL]|20UY[ l, and 2 = 3.7 | |Lee et aT.]|2Ull[ l. For details on 
how tne observed regions are determined, see tne main text. 



ing galaxies, the "SFR sequence", at different redshifts. 
This is related to the specific s tar formation rate ver- 
sus redshift discussed in Section [4.2[ but measured for a 
range of stellar masses at fixed redshifts rather than for 
a fixed stellar mass at a range of redshifts. 

To extract this quantity from our model grid, we pick 
a redshift slice to correspond to a chosen observational 
study, then plot stellar mass versus specific star forma- 
tion rate for all halos at that redshift. Figure [T3| shows 
the result, compared to several observational surveys: a 
sampl e of L yman Break Galaxies at z w 3.7 from |Lee] 



et al. pOll |), a sample of BzK galaxies at z w 2 from 
Daddi et al. (2007), a sample of 2 « 1 galaxies from 
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the GOODS survey studied by [Elbaz et~aL] ( |2007| ) , and 
z < 0.1 galaxies from the Sloan D igital Sky Survey stud- 



ied by Brinchmann et all (2004 ). For the Daddi et al 
( |2007| ) and |Lee et al.| ( |2011f samples we use the ob 



servational constraint region s giv en in |Lee et al.| fl2011|) . 
For the |Brinchmann et aL] ( |2004[ ) and [Elbaz et al.| (|2007|) 
sampl es, we use the fitting formulae given in |Elbaz et aT" 
(2007); since they do not give explicit confidence regions 



or stellar mass ranges, we set our mass range based on 
eyeball estimates of the region occupied by the data. 
Similarly, w e adopt a scatter of 0.4 dex for the Elbaz 



et al. 



( 2007) sample and 0.5 dex for the Brinchmann et a 



] 2U04 ) sample, based on an eyeball estimate of the scat 
term the data. We correct all observations to a lChabrierl 
fl2005|) IMF, which is a factor of 1.59 lower in SFR than 
• Kennicuttl ( fl998| > IMF. 



The plot shows that our fiducial model does a rea- 
sonably good job of reproducing the observed sSFR se- 
quence, particularly at high redshift. At lower redshift, 
our models fall slightly below the observed range, consis- 
tent with the simi lar u ndershoot of sSFR we saw below 
z = 2 in Section |4.2| Nonetheless, the level of agree- 



ment is gratifying given the extreme simplicity of our 
model. In particular, we recover the basic results that 
the sSFR is roughly independent of stellar mass at all 
redshifts, and that the value of the sSFR varies fairly 
slowly with redshift above z ~ 2. The model also makes 
clear predictions for the slope and normalization of the 
SFR sequence at higher redshift than has yet been ob- 
served. We do not show the corresponding predictions 
for our variant model with f g ,i n = 1, but, not surpris- 
ingly based on the discussion in Section [4.2[ it provides a 
much worse fit to the higher redshift data (and a some- 
what better fit to the low redshift data). The failure at 
high redshift is because of the lack of stellar accretion, 
which produces systematically smaller stellar masses at 
the same star formation rate. 

4.4. The Mass-Metallicity Relation 

Another observational constraint against which we can 
compare our model is the mass-metallicity relation. This 
is not tremendously sensitive to the inclusion or exclusion 
of metallicity-dependent star formation, because most of 
the galaxies for which this relation has been observed are 
in the regime where there is little mctallicity quenching. 
However, it is important as a consistency check. Since 
the model depends on metallicity affecting the SFR, it 
is critical that the metallicity evolution the model gen- 
erates be at least roughly consistent with the observed 
metallicities of galaxies as a function of mass and red- 
shift. 

Such an observational comparison is unfortunately 
very difficult because of uncertainties in the absolute 
zero points of the strong-line mctallicity indicators that 
are generally used to infer galaxy metallicities. While 
the various metallicity measurements in use can be cali- 
brated against each other to a scatter of ~ 0.03 dex, the 



zero point is un certain by as much as 0.7 dex (Kewley 
& Ellison 2008). Since our models predict an absolute 
metallicity, i.e. a mass of metals divided by a mass of 
gas, this means that there will be a large systematic un- 
certainty in any observational comparison. We handle 
this problem by fixing the models to match observations 
at a particular stellar mass and redshift, thereby fixing 




Fig. 14. — Mass-metallicity correlation for our fiducial model 
(solid lines) at z = 0.07 (green), z = 2.26 (red), z = 3.7 (blue), 
and 2 = 6 (purple), from top to bottom. The z = 3.7 and 2 = 6 
lines end at the highest stellar mass contained in our model grid 
at that redshift. For compa rison we show observed masses and 
metallicities: a sample from IT remonti et al.1 (12 004b at z = 0.07 
(green squares), a sample from hj'b et al. ( 21)06} at z = 2.26 (red 
circles), and a s ample from jM aiolino et al. (2008) at z = 3.7 (blue 
triangles). The |Tremonti et al. (2004J points are averages over 
many galaxies in tne S1JSB; tne error bars in mass indicate the 
mass range sampled by each data point, and the error bars in 
metallicity indica te the range from t he 16th to the 84th percentile. 
The points from |Erb et al.l (12006b are a sample of UV-selected 
galaxies; the error bars m mass indicate the standard deviation 
of the masses contributing to that point, and the error bars in 
metallicity indicate the errors in the inferred metallicities. The 
points from |Maiolino et al.| | |2008^ are individual galaxies, with the 
error bars in mass and metahcity indicating the uncertainties in 
those quantites for each galaxy. The observations have all been 
corrected to the same IMF and metallicity calibration scheme, and 
the conversion between metallicity and 12 + log(0/H) from the 
models has been fixed empirically. See the main text for details. It 
is important to note that, because of the systematic uncertainties in 
the strong line calibrations, both the slopes of the various observed 
relations and the offsets between them should be regarded as highly 
uncertain. 

the zero point. We can then compare the mass and red- 
shift dependence of the models to the observations. It is 
important to note that we are not looking for more than 
very rough consistency here, because the slopes of metal- 
licity versus mass and metallicity versus redshift depend 
on which metallicity calibration one adopts. Different 
choices of calibration can make a significant difference to 
how well the data fit the model. 

Figure [14] shows a comparison between the mass- 
metallicity relatio n in our models a nd a compilation of 
observations from Tremonti et al.| J2004|) at z = 0.07, 



Erb et al | (|2006|) at z 

at z = 3. 



2.27, and |MaioImo et al.| ( |2008 l 
"Xs with other comparisons, we generate 
the model predictions by finding the redshift slice in our 
model grid nearest to the target value and plotting mass 
versus metallicity at that sl i ce. W e correct all the ob- 
served masses to a Chabrier ((2005) IMF, and all the ob- 
served metallicitie s to the KpQ2 calibration scheme using 
the conversions of |Kewley &: Ellison (2008). We fix the 
absolute calibration of the models~5y forcing agree ment 
between the models and the Tremonti et al. 
at z = 0.07 and logM M1 = -0.5 
adopting a conversion 



_ ([2004] data 
This corresponds to 



12 + log(0/H) - 10.75 + \og(M z /M g ), 



(30) 

where Mz and M g are the masses of metals and gas, 
respectively. 
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Fig. 15. — Stellar mass versus H2 to stellar mass ratio for our 
fiducial model (solid lines) at 2 = (green), 2 = 1 (lavender), 
z = 2 (red), and 2 = 4 (blue). For com parison, at 2 = we show 
the results of the COLD GASS project | |Saintonge et al.|2011a|b[ |. 
Filled circles represent detections, and downward arrows indicate 
5cr upper limits from non-detections. The high-redshift predictions 
are to be tested with ALMA. 

As the Figure shows, the agreement between the ob- 
served and model metallicity evolution is generally good 
at high z. In the present-day universe our model pro- 
duces a somewhat steeper mass-metallicity relationship 
than the observed SDSS sample. While this could indi- 
cate a real defect in the models, it could equally well be 
a result of the particular metallicity calibration scheme 
we have used, since the slope of the SDSS sample (and 
the other samples) can vary by factors of ~ 2 depending 
on the calibration scheme used. If our models really are 
too fiat compared to reality, this may be related to our 
models' tendency to somewhat underpredict the star for- 
mation rate in sma ll gal axie s at low redshift, as we will 
see in Sections 4.2 and 4.3 However, since the agree- 
ment between our models and the observations is good 
at z > 2, our results should be robust there. 

4.5. The Hi Content of Galaxies 

Another feature of our model is its ability to predict 
the H 2 content of galaxies as a function of redshift and 
other galactic properties. Figure [15] shows the H2 mass 
fraction versus stellar mass predicted by our models. For 
comparison, at z = we s how the results of COLD GASS 
( Saintonge et al. 2011a|b ), a volume-limited survey of CO 
line emission (used as a proxy to infer H2) in several hun- 
dred galaxies at distances of 100 — 200 Mpc, selected to 
have stellar masses A/* > 10 10 M . As the plot shows, 
we find generally good agreement between our model and 
the observed amount of H2 in galaxies at z = 0, albeit 
with quite a large spread in the data. We recover the 
trend that the H2 fraction is (very weakly) declining with 
stellar mass. We also note that the data contain a signif- 
icant number of non-detections. Some of these represent 
star-forming galaxies with relatively modest molecular 
content whose true value of Mh 3 /-M* is likely only a bit 
below the reported upper limit. However, some also rep- 
resent quenched galaxies that are almost entirely devoid 
of star formation and molecular gas. Our model does 
not correctly produce galaxies of this type due to our 
assumption that accretion cuts off sharply at a single, 
relatively high, halo mass. The observed co-existence of 
star-forming and passive galaxies at the same stellar (and 



presumably halo) mass indicates that this is clearly an 
oversimplification. 

At higher redshift, we predict that the H2 fraction at 
fixed stellar mass should rise, reaching ~ 10% by mass 
at z = 1 and roughly 50% by mass at z = 2. This is 
qualitatively consistent with the observed trend toward 
high molecular gas fractions at z ps 2 (e.g. Tacconi et al. 



2010 1, but given the selection biases inherent in the cur 
rent high z data, it is hard to draw any general con- 
clusions yet. However, these predictions will be testable 
with ALMA in the next few years (e.g. see reviews by 
Combes|M0||2nTTl ). 



4.6. Stellar, H I, and H2 Mass Functions 

Another useful quantity to plot from our models are 
the mass functions of stars, H 1, and H2 as a function 
of epoch. Such a comparison can provide some physi- 
cal insight into the behavior of the model. However, we 
strongly caution that we should not expect good agree- 
ment between the model and observations here. Our 
models are very simple. Unlike in a full semi-analytic 
model, we have not made any effort to tune parameters to 
reproduce observations, and we omit much physics that is 
known to be necessary to reproduce observed mass func- 
tions. For example, we assume that there is one galaxy 
per halo at all halo masses, clearly not a reasonable as- 
sumption for high mass halos that should host clusters. 
We make the mass loading factor of galactic winds inde- 
pendent of halo mass or other properties, which is prob- 
ably not reasonable for low masses. Finally, we assume 
that all halos at a given mass and redshift are the same, 
ignoring scatter, which is not reasonable given the steep- 
ness of the halo mass function. 

It is important to note that these omissions are un- 
likely to dramatically affect the other observational com- 
parisons we present. In particular, we are not particu- 
larly concerned with galaxies larger than L*, and these 
do not contribute greatly to the cosmic star formation 
rate budget, so our incorrect treatment of them is not 
that important. Similarly, our neglect of scatter can 
have dramatic effects on the mass function, but it is un- 
likely to affect correlations between, for example, mass 
and metallicity. The one factor that could have a signif- 
icant effect is our assumption of weak galactic winds at 
low masses, since these can have much the same effect 
as metallicity quenching in suppressing star formation in 
small galaxies. In reality both metallicity quenching and 
strong winds are likely to be important, but we have de- 
liberately kept the winds weak and mass-independent in 
order to perform a cleaner experiment on the effects of 
metallicity alone. 

To construct our mass functi ons, we compute the halo 
abundance n(Mh,z) from the Sheth fc Tormen (19991 
approximation as in Section |4.1[ and for each halo we 
assign the Hi, H 2 , and stellar mass of the corresponding 
halo in our grid. Since the stellar mass is a monotonic 
function of the halo mass in our models, at any redshift 
Z the number density of galaxies n(M*,z) in the mass 
range from Af* to + dM* is simply given by 



n(M*,z) = n(M h ,z) 



dM* 
dM h 



(31) 



where on the right hand side Mh is the halo mass for 
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Fig. 16. — Stellar, H I, and H2 mass functions (left to right columns) at redshifts z = 0, z = 0.9, and z = 2.5 (top to bottom rows). In 
each panel, the thick black line is the fiducial model. The thin black line is the mass function that would be produced if all halos had a 
stellar, H I, or H2 mass equal to the halo mass times the universal baryon fraction /;,, e.g. in the left panels the thin black lines show the 
outcome of baryons accreting onto halos and converting into stars with perfect efficiency, in the middle panels the lines show the result of 
baryons accreting with perfect efficiency an d remaining enti rely H I, etc. Gra y lines and symb ols indicate observat ions. The stellar mass 
function points are from 



|Borch et ah] ( |2006| z = 11.8 — 1; crosses T TIPert 
z = 2 — 3; stars J, and [Mar chcsini ct al. (20(jy 



Li & White| ( |2009|l, as corrected by|Guo et aI. N2010| z 
TT crosses ), |Perez-Gonzalez et al.] 



( 2UUS 



"0.8 



- 0; asterisks), |Bundy et al.| j2006| 2: = 0.75 - 1; c i rcles), 
1, 2 - 2.5, and 'Z.b - 3; piusses), Fontana et al. (2006 
1 2UUb'| z = 2 — 3; triangles). All stellar data h ave been corrected to a |Ch.aprier| ( |2UUb| l IMF, and 
error bars have been omitted tor clarity, i'he z = H I mass function is the fit given by |Martin et al.| |2010| . Panels witnout gray lines or 
symbols indicate that no data is available. 



which the stellar mass is M*, and the derivative is eval- 
uated as this mass. 

The atomic and molecular gas masses are not strictly 
monotonic in the halo mass in our model, so there may 
be multiple halo masses M/, 1, 2l ■ ■ ■ for which the 
corresponding atomic or molecular mass is Mm or Mh 2 . 
The equivalent expression for the number densities of 
galaxies with a given H 1 mass is therefore 



n{M m ,z) =^2n(M h]i ,z) 



dM h 



(32) 



where the sum runs over all halo masses M^^ for which 
the H 1 mass is Mhi- The expression for H2 is analogous. 

Figure 16 shows the model mass functions at three red- 
shifts, compared to data where it is available. The plots 
reveal several interesting features of the model. First ex- 
amining the stellar mass function at z = 0, we see the 
sharp drop at low stellar masses corresponding to halos 
where star formation begins to be suppressed by metallic- 
ity effects. The mass at which this cutoff occurs moves to 
smaller values at higher redshifts. Well above this cutoff, 
our model provides very little suppression of star forma- 
tion, which is not surprising given that we have not in- 
cluded strong galactic winds or similar mechanisms, and 
that we do not have any metallicity scatter in out mod- 
els. Thus we suppress star formation too much below 
the cutoff mass and too little above it. This shows up in 



reverse in the H 1 mass function, where we overpredict 
the number of low H 1 mass galaxies, and underpredict 
the number of large H 1 mass galaxies. Clearly a sharp 
cutoff in star formation such as the one produced by our 
model is at best a crude approximation to reality, and 
the real effects of metallicity should be spread out much 
more in halo mass. 

We also fail to produce galaxies with very large H 1 
masses, indicating that we are probably converting H 1 
to H 2 and thence into stars too efficiently in some cases. 
The origin of this failure is not entirely clear. One pos- 
sibility is that it comes from our lack of scatter in spin 
parameter. Galaxies with large H 1 masses tend to have 
very large, extended H 1 disks that are non-star-forming. 
Our model does not produce these. 

4.7. The Cosmic Star Formation History 

By combining the star formation rate as a function of 
halo mass and redshift from our model grid with a cal- 
culation of the halo abundance as a function of mass and 
redshift, we are able to produce the cosmic star forma- 
tion history for our toy model. We caution that, since 
the agreement between our simple model and observed 
sSFRs at redshifts < 2 is only approximate, we should 
expect no better agreement here. Our goal is simply to 
understand the qualitative effect of metallicity quench- 
ing on cosmic star formation history with an emphasis 
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Fig. 17. — The contribution of a halo at a given mass and redshift to the cosmic star formation rate or accretion rate density, 
form; Mg ,i n ] {Mh , z^n^M^ , z)A4^ , indicated by color. The format is similar to Figure |?J Heavy black contours 
show the loci of halos and redshifts with dp^ ^g^i/dlnM^ = —2.5, —2.0, —1.5, and —1.0 Mq yr^ 1 Mpc -3 , as indicated. 



on the high redshifts where its effects are greatest. 

We c omp ute the halo number density n(A// l ,z) as in 
Section |4.1| The cosmic star formation rate density at 
redshift z is then given by 



p* 



M, .form 

(M h ,z)n(M h ,z)M h dlnM h . (33) 



The cosmic density of gas accretion into galaxies is given 
by the analogous expression with M*j orm replaced with 



M„ 



Note that this expression omits stars formed in 



the accretion shock when / 9 ,i n 1, but that is a modest 
effect. 



We plot the integrand of equation ( 33 ) , and the analo 



gous accr etio n rate for the cosmic gas accretion density, 
in Figure [17} The figure illustrates several points. First, 
our fiducial model range of masses does a good job of 
sampling the halo masses that dominate the star forma- 
tion rate, and thus we can make an accurate estimate of 
p* from it for the fiducial model. Second, it is obvious 
upon inspection that, for the fiducial model, the star for- 



mation rate density of the universe will peak at z ~ 2 — 3 
and fall off sharply on either side of it. On the low z side 
this fall-off is driven by the decline in star formation rate 
(which is in turn driven by the falloff in gas accretion 
rate). On the high z side, the decline in star formation 
rate is driven by the suppression of star formation due to 
the difficulty in forming a cold, molecular phase in small 
halos. In comparison, if we ignore the need to form such 
a phase and set /h 2 = 1, or if we were to assume that 
halos could form stars as quickly as they could accrete 
gas, the cosmic star formation history would be much 
less peaked toward z ~ 2 — 3. At high redshift it would 
be dominated by small halos undergoing rapid accretion 
and star for mation. 

Figure [l8| shows the cosmic star formation rate density 
as a function of redshift integrated over all halos, both 
as actually observed and in our variant models, and the 
cosmic accretion density into halos above a given mass. 
These latter lines show what we would expect if halos 
simply turned gas into stars as quickly as they accreted 
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Fig. 18. — Cosmic density of star formation in our fiducial model 
(solid red line), in the model where we set /jj 2 = 1 (dashed red 
line), in the (j = 0.8 model, and total gas inflow rates into halos 
with masses above logMj, yi = — 3, —2, and —1 (black lines, from 



top to bottom). We do not show the f g 



1, M ret = 0.03, or 



Esf = 10 models, because they differ only slightly from the fiducial 
model, and we do not show the Zigm = 2 x 10 — 4 model because 
it is qualitatively similar to the (j Q = 0.8 one. We also show an 
estimate of the observed cosmic star formation rate density (blue 
shaded region). The region plotted is a fit based on combinin g 
the SFR density es t imates compiled by Hopkins & Bcacom J2006L 
|Kistler et aT] l |2009[ l, |Bouwens et a,l.| ( |20lU] ) , and |Horiuchi fc ijeacom| 

it from the IGM; the model with a minimum halo mass 
Mh,i2 = .1 is a rough appro ximation of the model pre- 
sented by ( Bouche et al. 2010[ ), although it is simplified in 
that it assumes a star formation time of zero as opposed 
to som e number of galactic orbital periods as |Bouche| 
let al.l assume. 

As the plot shows, all the models recover the observed 
star formation rate from z = 0—2 reasonably well. Given 
the long timescales from z = 2 to the present, the star 
formation rate in galaxies at this epoch becomes limited 
only by their gas supply, and it is the decline in gas 
supply from the z = 2 to the present that drives this 
trend. Note that our model SFR is slightly above the best 
fit observations at z = 2, but the exact shape and value 
of the peak in our models depends on how we choose to 
truncate cold gas accretion in massive halos near z — 2. 

In contrast, the only models that are able to reproduce 
the observed star formation rate at z = 2 — 8 are those 
that include metallicity quenching (fiducial, f g> i n — 1 and 
Mret = 0.03) and the ad-hoc model with a thresho ld for 
accretion at Mh,vi, — 0.1 (as in Bouche et al.|2010 ). The 
run with /h 2 = 1 provides a more gradual rise of the SFR 
density at high redshift, similar to the ad- hoc model with 
a threshold for accretion at Mh,i2 = 0.01, indicating that 
the suppression because i acc is shorter than tsF is effec- 
tive below a characteristic halo mass of ~ 1O 1O M0. This 
is indeed comparable to the Press-Schechter typical halo 
mass at the end of this period, z ~ 2, where the SFR fi- 
nally catches up with the IR|jIn our fiducial model and 
its variants that also include metallicity quenching, we 
recover the steeper rise in time because star formation at 
high redshift is also suppressed by the difficulty in form- 
ing a cold molecular phase in small galaxies that have 

7 In fact, the model with /h 2 = 1 should have an even higher star 
formation rate density at high z than what is plotted, because our 
model grid does not fully sample the small halos that dominate the 
star form atio n rate at high z if we set /jj 2 = 1. This is illustrated 
in Figure 1171 



not yet enriched themselves with metals. The £i = 0.8 
model is in between the fiducial model and the /h 2 = 1 
model, but is still marginally consistent with the obser- 
vational constraints. The Zigm = 2 x 10~ 4 model (not 
shown in the plot) is very similar. Thus we see that the 
amount by which star formation is suppressed when we 
adopt a realistic treatment of ISM thermodynamics is 
somewhat but not tremendously sensitive to how metals 
are retained in halos as a function of mass and to the 
metal enrichment history of the IGM. 

In the simplified model where the star formation rate 
is set equal to the accretion rate into halos with masses 
below Mh,\2 = 0.1, the mass limit has much the same 
effect as including metallicity in the star formation law 
on the integrated star formation rate of the universe. 
Both suppress star formation at high z. This demon- 
strates that quenching below a threshold mass captures 
some of the effects of metallicity-dependent star forma- 
tion. The model presented here provides a physical mech- 
anism that induces a strong but more gradual mass de- 
pendence, which boils down to a similar effect on the 
cosmic SFR history. Our physical model, while it does 
reduce the SFR in small halos, does not make them en- 
tirely devoid of stars, and it still allows galaxies like the 
Small Magella nic Cloud which lives in a h alo of mass 
M h ,n ~ 0.01 (Be kki fc Stanimirovic]|2009|) . In reality, 
the suppression ot star tormation is probably even less 
sharp than in our models, since there are almost certainly 
large stochastic variations in the metal enrichment his- 
tory of small halos that our simple deterministic model 
does not capture. The difference between our fiducial 
model and the (i a = 0.8 variant suggests what effect this 
stochastic variation is likely to have. We should therefore 
emphasize that the value of the effective threshold mass 
is predicted by our model only as an order-of-magnitude 
estimate. 

5. CONCLUSION 

In this paper we investigate how a metallicity- 
dependent star formation rate affects the evolution of 
star formation in growing galaxies, and how it impacts 
the observable cosmic star formation history. The physi- 
cal basis for the effect is that star formation occurs only 
in the molecular phase of the ISM where the gas is able 
to reach very low temperatures and low Jeans masses. 
Because metal atoms and dust grains are the primary 
coolants and catalysts of H2 formation, and because dust 
grains are also the primary shield against both photo- 
electric heating and H2 dissociation by UV photons, the 
ability of the ISM to form a cold molecular phase depends 
strongly on its metallicity. 

In order to study how these effects modify the star for- 
mation histories of galaxies, we have developed an ideal- 
ized model in which we followed the evolution of the mass 
and surface density of a galaxy as it accretes baryons, 
forms stars, and self-enriches with metals. The model is 
not meant to provide a perfect fit to observations, since 
we have omitted numerous important processes, but it 
does allow us to determine qualitatively how the phase 
transition from warm H I to cold H2 modifies the observ- 
able features of cosmological star formation history. The 
basic physical model ingredients are: 

1. Dark matter and baryons flow into galaxies in an 
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average cosmological rate including mergers. 

2. Disks have exponential surface density profiles, 
with scale lengths determined by cosmological evo- 
lution under conservation of angular momentum 
and the assumption that dark matter and baryons 
share the same spin parameter. 

3. Mass is conserved in galaxies as baryons flow in, 
gas accumulates, turns into stars, and is lost to 
outflows. 

4. The star formation l aw is that proposed by 



Krumholz et al. 
area is close to 



((2009cJ), in which the SFR per unit 
inearly proportional to the column 



density of the molecular phase of the ISM. 

5. The characteristic surface density at which the gas 
switches to the cold, molecular phase is roughly 
given by E ~ IO/(Z/Z )- 1 M pc~ 2 , and near 
this threshold the cold molecular fraction roughly 
obeys /h 2 cx S(Z/Z ). 

6. Metal enrichment is a function of the SFR, and the 
fraction of metals retained in a galaxy rather than 
lost to the IGM is a function of halo mass. 

7. The fraction of ex situ stars in the accretion flow is 
determined self-consistently over the cosmological 
halo population. 

We find that at very high redshifts, z > 4, the baryonic 
mass of the main-progenitor galaxy is growing rapidly 
with time (roughly cx e~ z ), with an inflow rate that at 
z ~ 8 is already almost as high as at z ~ 2. However, 
at these early times the SFR cannot catch up, because 
the timescale for star formation in the dominant mode 
is longer than the timescale for inflow. Furthermore, the 
low metallicity due to the short history of star formation 
and the small escape velocity in the typical high- z galax- 
ies is responsible for a low molecular fraction, and this 
further lowers the SFR compared to the inflow rate. The 
resultant SFR is growing steeply as exp(— 0.65z) until 
z ~ 2, contrary to the popular model assuming an ex- 
ponential decay of SFR. Most of the accreted gas during 
this period accumulated in an atomic gas reservoir, wait- 
ing for the metallicity to grow before it could produce a 
cold molecular phase and turn into stars. 

As time passed in the models, the ratio of timescales 
for star formation and inflow declined quite rapidly, the 
forming stars gradually enriched the ISM with metals, 
and the growing halo mass became more capable of re- 
taining them. The SFR thus grew rapidly, and caught 
up with the IR by z ~ 2. This is cons i stent with ob- 

When 



(2010). 



servational findings by |Papovich et al. 
integrated over the halo population at a given time, the 
cosmic SFR density history (SFH) rose with time in the 
models, in agreement with the observed trend. The in- 
clusion of metallicity in the star formation process affects 
the SFH in a way that is similar to the effect of a sharp 
threshold for star formation at halo mass M t h ~ I0 11 M . 

We caution that the exact value of the halo mass at 
which metallicity quenching begins to suppress star for- 
mation is not well-determined in our models. It depends 
on the overall life cycle of metals in galaxies and the IGM, 
and is sensitive to quantities such as the metal yield of 



supernovae, the ability of small galaxies to retain the 
metals they produce, and the metallicity of the IGM gas 
that accretes onto galaxies at high z. Thus our conclu- 
sions about how metallicity quenching alters the cosmic 
star formation history should be regarded as qualitative 
rather than quantitative. However, the qualitative ef- 
fect that star formation is suppressed below a threshold 
mass is robust, and appears in all our models regardless 
of how we alter the physics of the metal cycle. This has 
the effect of suppressing star formation in small galaxies 
at high redshift, and shifting some of this star formation 
to lower redshift after metals have had a chance to build 
up. The exact amount of suppression and shifting de- 
pends on the exact threshold for metallicity qu enching. 
A threshold of this sort was indeed advocated by|Bouchi 



et al. ( 2010| their Fig. 1) as a possible ad- hoc explanation 
for the SFH rise, as well as several other phenomena in- 



cluding the downsizing of galaxies ( Neistein et al.| 2006 ) 
The suppressed in situ SFR makes the stellar fraction 
in the galaxy small, and the growth of stellar mass dom- 
inated by accretion of stars that formed ex situ or in 
mergers. Under this condition, we showed that the sSFR 
tends to be roughly constant in time, from very early 
times till z ~ 2. This is consistent with the observa- 
tional indications for a sSFR plateau for g alaxies of a 
fixed mass o bserved in the range z = 2 — 8. IWeinmann 



et al. (2011) showed that this plateau is rather puzzling 
when compared with the steeply declining inflow rate 
with time, but this work shows that including the metal- 
licity dependence helps resolve the issue. One caveat in 
our otherwise self-consistent model is that we were forced 
to select a small initial stellar fraction for halos at early 
times, trying to mimic starbursts in low-metallicity high- 
density conditions induced by mergers. Although the ex- 
istence of a sSFR plateau is not sensitive to this fraction 
(as long as it is non-zero) , the quantitative value and du- 
ration of the sSFR plateau is, as described in equation 
27j sSFR - 2(/*a n /0.I)- 1 Gyr" 1 , where /*, in < I is the 
Traction of the inflowing baryons at high z that either ar- 
rive already as stars, or that turn into stars in a starburst 
immediately upon accreting. 

We note that the dominance of ex situ over in situ SFR 
at high redshift is the opposite of what was suggested 



by the simulations of Oser et al. (2010). While their 
simulations had efficient SFR in small galaxies at high 
redshifts, the inclusion of the effects of metallicity on 
the SFR introduces strong suppression of in situ star 
formation in halos < 10 11 M at z > 2. This allows a 
small fraction of ex situ stars to dominate, and naturally 
lead to a sSFR plateau. 

The gas reservoir that builds up at z > 2 because it 
cannot form a cold, molecular phase capable of turning 
into stars provides additional fuel for star formation in 
10 10 - 10 12 M Q halos at z ~ 1 - 3, when the metallic- 
ity reaches values high enough to convert the bulk of the 
ISM to molecules. This allows the SFR to exceed the in- 
stantaneous gas IR, which helps to explain the otherwise 
puzzlingly high SFR density and sSFR observed in this 
redshift range. The reservoir gas can also help explain 
the observed massi ve outflows from z ~ 2 — 3 galaxies 
tiSteidel et al.|2010[ ), which may exceed the inflow rate of 
fresh gas minus the SFR. 

We specify additional observable predictions of our 
model, as a function of galaxy stellar mass and redshift. 
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First, the sSFR sequence, of sSFR versus stellar mass, is 
rather fiat below M* ~ 1O 115 M at all redshifts. This is 
complementary to the constancy of the sSFR with time 
at z > 2. Second, at any redshift, the metallicity is ris- 
ing roughly oc m} /2 in the range M* < 1O 1O - 5 M , and 
it flattens off at higher masses. At a given stellar mass, 
in the range z > 2, the metallicity grows in time roughly 
oc z~ 3 . Third, the H2 fraction is very weakly decreasing 
with mass at fixed redshift, and is higher at higher red- 
shifts. At z w 2, H 2 masses are roughly equal to galactic 
stellar masses. We should add that the slow SFR and 
the resultant gas accumulation at high redshift is likely 
to lead to extended discs with low bulge-to-disk ratios, 
and thus help explain one of the most interesting open 
questions in galaxy formation ( Guedes et al.|20iT Brook 
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APPENDIX 

A. PARAMETERS FOR GAS EXPULSION AND STELLAR RECYCLING AND YIELD 



We adopt e ou t = 1-0, based on the estimates of Erb (2008) that galactic winds tend to carry mass fluxes comparable 
to the galactic star formation rate. In reality the wind mass loading factor almost certainly depends on galactic 
properties such as the halo circular velocity, star formation rate, or star formation rate surface density. We do not 
attempt to include this dependence. In part this is because it is extremely poorly determined both observationally 
and theoretically. However, experimenting with different values of e ut, including mass-dependent ones, shows that 
is has very little effect on metallicity quenching. Halos that are small enough to have their star formation quenched 
by metallicity effects do not form many stars,, and thus it matters fairly little what their mass loading factor is. 
Conversely, since galactic winds that do not preferentially carry away metals (i.e. those described by e out rather than 
C in our formalism) do not alter the metallicity, the choice of e ou t changes the absolute stellar mass of large halos, but 
it does not change the extent to which t hey are affected by metallicity quenching. 

For R and y, following the notation of iTinsleyl (fl980|), we define the return fraction by 



R = 



f™" [m — w{m)\ 4>{rn) dm 



dm 



(Al) 



where <j)(m) is the stellar IMF, w(m) is the mass of the remnant left by a star of mass to, and mi and m u are the lower 
and upper mass limits of the IMF. The yield of element i is defined by 



1 Imi mpi(m)(f)(m) dm 
1 — R f™ u m<p(m) dm 



(A2) 



where Pi(m) is the net fraction of the star's initial mass that is converted to element i in a star of mass m. 
We adopt the instantaneous recycling approximation, and set 

r to, m/M < 1 

w(m) = \ 0.7 M Q , 1 < to/M© < 8 , 
( 1.4 Mq, m/M Q > 8 



(A3) 



corresponding to an assumption that stars of mass m < 1.0 Mq never leave the main sequence, those of mass m = 1 — 8 
M & pr oduce white dwa rfs of mass 0.7 Mq, and those of mass m > 8 Mq produce 1.4 Mq neutron stars. For </>(m) we 
use the Chabrier (2005) form, 



(to) = 



0.41m" 
0.18m" 



exp 
2.35 



(logm-log0.2) J 



2x(0.55) 2 



TO < 1 Mq 
TO > 1 Mq 



(A4) 



within the range to./ — 0.08 Mq, m u = 120 Mq. With this IM F and value of w(m), we have R = 0.46. 

For the production function Pi(m), we adopt the models of Maederj ([1992) for Z = 0.001 (i.e. low metallicity) stars 
(his Table 5). Since we are treating metallicity with a single parameter, we use the total metal yields, and interpolate 
between the tabulated values. With these models and our chosen IMF, the yield is y — 0.069. If we instead use the 
models for Solar metallicity, y = 0.054 instead, and the change in results is small. We prefer to use the low-metallicity 
value because our models are most sensitive to the metallicity evolution when the overall metallicity is small. 

One might worry that instantaneous recycling becomes a poor approximation in the early universe, where stellar 
lifetimes are not n ecessarily short compared to the Hubble time. H owever, the approxim ation is surprisingly good. 
Using starburst99 ( |Leitherer et al. 1999 Vazquez & Leitherer 2005) for a |Kroupa| (2002) IMF with an upper mass 
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cutoff of 120 Mq, we find that the return fraction reaches R = 0.27 after only 100 Myr, i.e. about 60% of the mass 
that our instantaneous recycling approximation predicts will eventually be returned to the ISM is returned in the first 
100 Myr after star formation. Metal production is even faster. For comparison, at z = 30, the redshift at which we 
begin our calculations, the age of the universe for our cosmological parameters is also very close to 100 Myr. Thus the 
error made by the instantaneous recycling approximation is at most a factor of ~ 2 very early in our calculations, and 
becomes significantly smaller than that even by redshift 10. 

B. COMPUTING THE GASEOUS INFLOW FRACTION 

The gaseous inflow fraction f g \ n is expected to depend on the halo mass Mh and redshift z. We estimate / 9i i n self- 
consistently using the grid of models we describe in Section [3] The procedure i s as follows. F i rst, w e must estimate 
how the total accretion rate is divided up among halos of different masses M- m . Neistein et al. (2006) (their equation 
A15) show that the total accretion rate onto a halo of mass Mh arising from halos of mass less than M; n may be 
approximated by 

M h = , Mh (Bl) 

V tt ^SpT) - S{M h ) 

where S{M) = a 2 {M) is the variance of the density fluctuation in t op- hat spheres contai ning average mass M; we 
approximate this usi ng th e fitting function given in equation (A16) of Neistein et al. (20061. The total accretion rate 
is given by Equation ( |B1[ ) evaluated with Afj n s» M/j/2, and the differential contribution to the accretion rate by halos 
in the mass range M- m to M- m + dM- m is given by 



dM h 
dM;„ 



dS(M in ) 



2tt [S(M- m ) - S(M h )} 3 / 2 dM ix 



(B2) 



Note that Mk is non-zero even f or M- ln = 0, because some accretion is contributed by matter that is not in an existing 
halo (see also Genel et al. 2010). We define the non- halo accretion rate 



M h (M in = 0). 



(B3) 



Numerical evaluation shows that Mh, n on-hsdo / Mh ~ 0.2, decreasing very weakly with Mh- Given these expressions, 
we can approximate / 9i i n for a given halo by computing the accretion-weighted mean gas fraction of all the halos that 
are accreting, i.e. 



9,m 



1 

~M h 



M h /2 



f a (M in )^dM in + Mh .non— halo 



(B4) 



Note that we have assumed here tha t the baryonic component of the non-halo accretion is pure gas. 
In practice, we evaluate Equation (B4) numerically as follows. We assume that there is no star formation in halos 



smaller than the smallest halo in our grid of 400 model halos (which is true at least for our fiducial model), and so 
we adopt f g (Mi n ) = / s ,mit for all values of M in smaller than the smallest halo in our model grid. Thus, the stellar 
fraction of these halos never changes. For larger values of M; n , we store our grid of models at 400 redshifts from 30 
(our starting redshift) to 0, and evaluate f g by interpolating on the grid in halo mass and redshift. This procedure does 
not require any iteration to converge, because for a halo of mass Mh at a given z, we only need to evaluate f g {M- m ) 
for values of M; n < Mh/2. Provided that we compute the evolution of our model grid starting with the smallest halo 
and proceeding upward in mass, by the time we reach a halo of a given mass, we already know the full evolutionary 
history of all lower mass halos, and thus those parts of the interpolation grid we require have already been filed in. 
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